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Foreword and introduction 

My collaboration with Gabriele Veneziano has continued, almost uninterrupt- 
edly, for more than fifteen years (even now we are preparing a joint contribu- 
tion to the book "Beyond the Big Bang", which will be published by Springer, 
as this book). Our first meeting dates back to 1989, when Gabriele came to 
the University of Turin to give a series of talks and seminars. At that time 
I was working there as a young researcher at the Department of Theoretical 
Physics, and I remember that Sergio Fubini (professor at the same Depart- 
ment) introduced me to Gabriele before a seminar. After the seminar we went 
to my office, and we started talking about cosmology, big bang, inflation, and 
strings. Gabriele was able to make me feel at ease, in spite of the fact that I 
was a bit embarrassed, being face to face with such a world-renown scientist 
like him: before that meeting, indeed, I knew him only for having seen his 
name quoted in many books and articles as one of the founders of string the- 
ory. I could not imagine that I was about to embark on the most stimulating 
and important adventure of my scientific life. 

After that meeting we started collaborating on string cosmology, and I 
visited very often the Theory Division (now "TH Unit" ) at CERN, living in 
Geneva also for long periods. This has given me the opportunity to appreciate 
Gabriele not only as a scientist - whose inventiveness, originality, profundity 
of thought will not be stressed here, because they are well-known to the physi- 
cist community - but also for his human qualities. His tireless enthusiasm for 
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physics, his generosity in sharing knowledge, his intellectual honesty, always 
make working with him a rewarding and enjoyable experience. I have count- 
less memories of days spent discussing and working out calculations on the 
blackboard of his office (see Fig. 1), with short "coffee breaks" every now and 
them, talking about physics even during lunch and dinner. Countless are the 
things I have learned from him, not only from a scientific but also from a 
human point of view. I will be grateful to him forever. 




Fig. 1. Gabriele Veneziano (left) and the author (right), talking about dilatons at 
CERN (January 1994). 



Choosing among the lines of research developed in collaboration with 
Gabriele, I will concentrate the contribution to this book on the possible 
role played by the dilaton in a cosmological context, with particular atten- 
tion to the phenomenological aspects of dilaton cosmology. The dilaton is a 
fundamental scalar field appearing in all models of superstrings, dilaton cos- 
mology is probably the most natural and typical form of "string cosmology" , 
and a direct/indirect confirmation (or disproof) of its predictions could give 
us important experimental information on string theory in general. 

The present contribution contains three lectures. The first lecture (Sect. 1) 
is devoted to the presentation of a primordial cosmological scenario in which 
the background evolution is dominated by the dilaton, and the Universe is 
driven through an accelerated phase representing the "dual" counterpart of 
the standard, decelerated evolution. The second lecture (Sect. 2) will discuss 
the possibility that a cosmic background of relic dilaton radiation could have 
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survived until present, and could be detectable by the gravitational antennas 
that are presently operating (or planned for a near- future operation) . Finally, 
the third lecture (Sect. 3) will suggest a possible "dilatonic" origin of the dark- 
energy fluid dominating the cosmic acceleration recently observed on large 
scales, stressing the main differences from other, more conventional models of 
scalar "quintessence" . 

Notations and conventions 

Unless otherwise stated, the following conventions are used throughout this 
paper: Greek indices run from to d, Latin indices from 1 to d, where d = D — 1 
is the number of spatial dimensions of the D-dimensional space-time manifold. 
The metric signature is: 

g^ = diag(+, •). 
The Riemann tensor and its contractions are defined by: 

r p — F) r p 4- p p r p — (n^iA 

The conventions for the covariant derivative are: 

v M u Q = d„v a - r m %, v M u« = d^v a + a vK 

Finally, we often use the convenient notation: 

(V0) 2 = V M 0V0, V 2 = V„V>. 



1 Dilaton-dominated inflation: the pre-big bang scenario 

If we apply to the "specialized" literature for a description of the birth and 
of the first moments of our Universe, we may read, in the (probably) most 
ancient and authoritative book, that 

"In the beginning God created the Heaven and the Earth, 

and the Earth was without form, and void; 

and the darkness was upon the face of the deep. 

And the Breath of God 

moved upon the face of the water. " 



(Genesis, The Holy Bible). 
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The most impressive aspect of these verses, for a modern cosmologist, is prob- 
ably the total absence of any reference to the extremely hot, kinetic, explosive 
state that one could expect at (or immediately after) the "big bang" deflagra- 
tion. The described state, instead, is somewhat quiet, dark, empty - we can 
read, indeed, about "void" , "darkness" , and "the deep" gives us the idea of 
something enormously desert and empty. In this static configuration there is 
at most some small fluctuation (the "Breath"), a ripple on the surface of this 
vacuum. 

It is amusing to note that a state of this type (flat, cold and vacuum, 
only ruffled by quantum fluctuations), can be obtained as the initial state of 
our Universe, in a string-cosmology context, under the hypothesis that the 
Universe evolves in a "self-dual" way with respect to the symmetries of the 
low-energy string effective action [I, 2]. 

To introduce this result we start considering the gravi-dilaton sector of the 
low-energy effective action. To lowest order in the a' (higher-derivative) and 
gl (higher-loop) expansion the action is the same for all models of superstrings 
[3, 4], and is given by 

S = -rir / d d+1 x^e^(R + d^d^ + V). (1) 

Here </> is the dilaton, and A s = (27ra') 1 / 2 is the fundamental length parameter 
of string theory. We have written the action using the so-called String frame 
(S-frame) metric, i.e. the metric to which a "test" string is minimally coupled 
and in which its evolution is geodesies. We have also included, for completeness 
and for further applications, a (possibly non-perturbative) dilaton potential, 

v = v{4>). 

This action should be completed by the source term S m (g,(j)), describing 
the matter- fields contributions, and by the Gibbons-Hawking boundary term 
Ss, which is required (as in general relativity) to cancel the variational con- 
tributions of the second derivatives of the metric following from the Einstein- 
Hilbcrt Lagrangian \J~ r g R. For the S-frame action (1) the boundary term 
takes the form [5] 

S * = ^3=1 / ^9 e~*K a dS a , (2) 

where K a = Kn a . Here K is the trace of the extrinsic curvature of the d 
dimensional hypersurface Oil bounding the hypervolume fl over which we are 
varying the action, and n a is the unit vector normal to this hypersurface. 

The variation of the total action S + S m + Ss with respect to g^ v leads 
then to the equations 

G„„ + V M V„0 + \g^ (V0) 2 - g^ 2 4> - \ 9llv V{4>) = Af- V T^, (3) 

where G^y is the Einstein tensor and the gravitational stress tensor of the 
matter sources, defined as usual by the functional differentiation of S m as: 
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The variation of the total action with respect to (j> leads to the dilaton equation 
of motion, 

f)V 

i? + 2V 2 0-(V0) 2 + T/-^ = A^ 1 e < V, (5) 

where a is the (S-frame) density of dilaton charge of the sources, defined by 
the functional differentiation of S m with respect to 4>: 

2 SS m 

Using the dilaton equation to eliminate the scalar curvature, present in the 
Einstein tensor, we can eventually rewrite Eq. (3) in the convenient (simpli- 
fied) form 



(7) 



1.1 Scale factor duality 



We will now consider the particular case in which the space-time manifold 
described by the S-frame metric is spatially flat, homogeneous (but not neces- 
sarily isotropic), and in which the matter fields can be phenomenologically de- 
scribed as perfect fluids, at rest in the comoving frame of the given Robertson- 
Walker geometry. We can thus set, in the synchronous gauge, 

= diag(l, -a 2 <%), a 4 = a,(t), cj> = cp(t), 
Tjx V = diag(p, —piSj), P = p(t), Pi=Pi(t), a = a(t).(8) 

Separating the time and space components of the gravitational equations we 
then obtain, from the (00) component of Eq. (3), 

<}> 2 -2^j2^+(j2 H ) -J2 H ?- v = 2X s~ 1 ^p (9) 

i \ i / i 

(where Hi = di/a$). From the (u) component of Eq. (7) we have: 
From the dilaton equation (5) we are lead, finally, to 

20-^+2^^-^(2^ + ^)-^^ +y-|K = As d - 1 cV (11) 
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We have thus obtained a system of d + 2 equations for the 2d + 3 unknowns 
{oi, 4>, p,Pi, cr}: its solution requires the input of d+1 "equations of state" , Pi = 
Pi(p), a = cr(p), specifying the properties of the considered matter sources. 

Let us now consider the symmetries of this system of equations. There 
are two symmetries, in particular, that are relevant for the discussion of this 
section. One of them (also present in the cosmological equations of general 
relativity) is the invariance under the time-reversal transformation t — > —t, 
which implies 

Hi^-Hi, i^i?,, 4>^4>- (12) 

Thanks to this invariance property, if the set of variables S — {<ii(i) , <j>{t) , p(t)} 
represents an exact solution of Eqs. (9)-(ll), then the time-reversed set 
S = {cii(—t),(f>(—t), p(—t)} also corresponds to an exact solution of the same 
equations (with different kinematic properties, in general). 

The string-cosmology equations, in the particular case a = and V = 
const, are also invariant under other transformations which have no analogue 
in general relativity, and which include the inversion of an arbitrary number 
of scale factors of the background geometry (8): the so-called "scale- factor 
duality" transformations [1, 6]. For a simple illustration of this property we 
may conveniently rewrite the equations in terms of the "shifted variables" 4>, 
p, p i , a, defined by 

cj) = 4> — hi]^[ en = <f> — lnaj, i = l,...,d, 

i i 

P = P~[[ai, Pk=Pk~[[ai, o : = crY[a t . (13) 

i i i 

Eqs. (9)-(ll) then become: 

f -Y / H?-V = 2\t 1 e*p, (14) 

i 

H i -Hj> + \ d ^=X^e* (15) 
2%~f -J2H? + V-^=\ d s - 1 e*a. (16) 

i 

Under the transformation a -^>a = a -1 , on the other hand, we have: 

-ida ~ __ida da -1 

H = a L — -> H = a — = a— — = -H. (17) 
dt dt dt v 1 

We can then easily check that Eqs. (14)-(16), in the particular case a — and 
dV/d<j) — 0, are invariant under the scale-factor duality transformations: 
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This type of transformation is called "dual" as it generalizes to the case of 
time-dependent backgrounds the T-duality transformation inverting the com- 
pactification radius (thus interchanging "winding" and "momentum" modes) 
in the spectrum of a closed string, quantized in the presence of compact spatial 
dimensions [7]. For the invariance under the transformations (18), however, 
there is no need of a compact geometry; what is required, instead, is a non- 
trivial transformation of the dilaton. Let us suppose, in fact, that we are 
inverting a number n of scale factors, say ai, . . . , a n , with 1 < n < d: the 
condition <f> — > <p then implies 

d d n d 

(/) — lna^ = (/) — lnSi = — lna,^ 1 — lnaj, (19) 

i— 1 i— 1 i— 1 i—n+1 

from which 

n 

4>^4> = (f>-2^2\iia l . (20) 

i=l 

In the presence of sources, their energy density is also non-trivially trans- 
formed: the condition p — > p implies, in fact, 

d n d 

P\{a l =p\{a- 1 J] a *> ( 21 ) 

i— 1 i—1 i— 1+n 

from which 

n 

P^P = p\{al (22) 
i=l 

The transformation of the pressure is similar, but with an additional "re- 
flection" of the equation of state along the spatial directions affected by the 
duality transformation: 

n 

Pi->Pi = -PiY[a 2 k , i=l,...,n. (23) 

fe=i 

In any case, given a set of variables S = {a,i(t), <f>(t), p(t),pi} representing an 
exact solution of Eqs. (14)-(16), a new solution can be obtained by inverting 
an arbitrary number n (between 1 and d) of scale factors, and is represented 
by 

S = {a^ 1 ,^ 1 , . . . ,a~ 1 ,a n+ i, ...,a d , <j>, p, Pi,--- ,Pn, Pn+i, ■ ■ ■ ,Pd}, (24) 

where (f>, p, pi are given by Eqs. (20), (22), (23), respectively. 

The invariance under the transformations (18) is only a particular case of a 
more general 0(d, d) symmetry of the tree- level string cosmology equations [8] 
(see also the contribution of Meissner [9] to this volume) , and can be extended 
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so as to include the NS-NS two-form B^ v in the effective action. Such an ex- 
tension is also possible in the presence of fluid sources: a homogeneous gas 
of strings, in particular, provides a realistic example of source which is auto- 
matically compatible with the 0(d, d) symmetry of the background equations 
[10]. 

In addition, the invariance under the transformations (18) can be extended 
to the case of non-trivial potentials, dV/d<p ^ 0, and non-zero dilaton cou- 
plings to the matter sources, a ^ 0. In both cases, however, we need to 
generalize those parts of the action describing the self-coupling of the dilaton 
and the dilaton couplings to the matter fields present in S m . 

In the case of the dilaton potential it is well known [8]- [10] that the in- 
variance under the transformations (18) holds for non trivial V , provided V 
depends on <f> through the variable <f). Such a variable, unlike 0, is not a scalar 
under general coordinate transformations (as evident from the definition (13)): 
it is thus impossible, in a generic background, to define a potential which is 
function of <j) and which can be directly inserted as a scalar into the covariant 
action (1). However, as first pointed out in [11], the action and the corre- 
sponding equations of motion can be written in a generalized form which is 
invariant under general coordinate transformations in any metric background, 
using for the potential a non-local variable which exactly reduces to <j> in the 
limit of a homogeneous geometry. 

Here it will be shown that the invariance under the duality transformations 
(18) can be restored also in the presence of the dilaton charge a, provided the 
dilaton coupling to the matter sources is parametrized by a non-local variable, 
as in the case of the potential. This result is new, and will be explicitly derived 
in the following subsection. 

1.2 Non-local dilaton interactions 

The formalism introduced in [11] is based on the non-local variable £(x), 
defined by 

£(z)=£[#c)]=-lny ^(vqe-*VW) 9 %-«, (25) 
where we have explicitly inserted the parameter 

•=*(ww={-l; ( ( v$< o; < 26 » 

so as to include in the formalism both time-like and space-like dilaton gradi- 
ents. Note that we are using the convenient notation in which an index ap- 
pended to round brackets, (. . .) x , means that all quantities inside the brackets 
are functions of the appended variable. Similarly, <p x = <p{x). We can immedi- 
ately check that, for a homogeneous background of the type (8) with spatial 
sections of finite comoving volume (J d d y = Vd = const < oo), the variable £ 
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exactly reduces to the variable 4> of Eq. (13). In that case, in fact, an explicit 
integration gives 

e = 0-InJJa i -In^V (27) 
and the constant volume factor can be simply absorbed by rescaling 0, so that 

Let us now suppose that the matter couplings and the self-coupling of the 
dilaton are both parametrized by £, according to the effective action 

S = - / d d+1 x V^o-f [R + (V0) 2 + V(e-t)] 

+ J d d+1 x J—g £ m (c-«) + S E , (28) 

which is a (generally-covariant) scalar functional of the non-local variable £. 
Note that, without loss of generality, we have written both the potential and 
the matter Lagrangian C m as a function of exp(— £). In higher-dimensional 
manifolds with compact spatial sections, in fact, the exponential of the shifted 
dilaton plays the role of a "dimensionally-reduced" coupling parameter, and 
we may thus expect (at least in a perturbative regime) that dilaton inter- 
actions appear as a power expansion (or as a a simple function) of such an 
exponential [11]. 

The generalized equations of motion can now be obtained by computing 
the functional derivative of the action (28) with respect to g^ v and 4>. The 
derivative with respect to the metric, using the standard definition of grav- 
itational stress tensor, Eq. (4), and the properties of the delta distribution, 
leads to the (integro-differential) equations of motion 

G>„ + V M V„tf> + hg^ (V<^ 2 - 2V 2 </> - V) 

- l^V^yWi^Iv - 2A s d - 1 / m ) = Af-V T„„, (29) 

which generalize Eq. (3) (see Appendix A for the details of the derivation) . 
Here 

TV = to* - (y0)2 . (3°) 



I v {x)=\- d J & d+1 y(J—gV') y S{ct> y -^), (31) 
I m (x) = \- d J d d+1 y(^0/(^-^), (32) 



where the prime denotes the derivative with respect to the argument cxp(— £), 
namely: 
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The functional derivative with respect to cf) leads to the dilaton equation of 
motion, 

R + 2V 2 - (V0) 2 + V + ^ V XJJ ^ - 2A ^ 1J ™) 
+ (c-« - e~*J) (V - 2\ d s - 1 c*C' m ) = 0, (34) 
generalizing Eq. (5). Here 

J{x) = K d J d d+1 y (V=gy/3yW) S'fa - 4> v ), (35) 

where 5' denotes the derivative of the delta function with respect to its argu- 
ment (see Appendix A) . The combination of Eqs. (29) and (34) finally leads 
to the equation 



+\g^ (e"« - e"* J) (V - 2A s d " V^) = A s d " V T^, (36) 

generalizing Eq. (7). 

We can easily check that these new equations, written for a homogeneous 
background, are invariant under scale-factor duality transformations even in 
the presence of non-trivial potentials and dilaton couplings, i.e. for dV/d(, ^ 0, 
dCm/d^ ^ 0. Consider, for instance, the background configuration of Eq. (8) 
with time- like dilaton gradients, for which e = 1. From Eq. (30) we obtain: 

7o° = 0, ll=Si. (37) 

The (0,0) component of Eq. (29) thus coincides with the (0,0) component of 
Eq. (3), and is given by Eq. (9), as before. 

For the spatial components we first note that, performing the homogeneous 
limit in which we are lead to the identities 

^W(c-*I v -2\t 1 I m ) = c-Z(V'-2\t 1 c'»£' m ) 

W" 2A ^) ! (38) 



E t H t fdV ^-l ^Crn 

,d(f> 3 d(p 



(39) 
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Using such identities we find that the dependence on V and C' m completely 
disappears from the spatial components of Eq. (36) with e = 1, and we obtain 
the condition 

i2i J ' + ViV J > = Af- 1 e*T < J. (40) 
Written explicitly, the new spatial equation is given by 



Hi — H t 



(t-Y^H^j =\t 1 e* Ph (41) 



and is thus crucially simplified with respect to the corresponding (local) spatial 
equation (10). 

The dilaton equation (34) also simplifies in the homogeneous limit, thanks 
to the identities (38) and (39) from which we obtain 

i? + 2V 2 ^-(V0) 2 + y-^= = -2Af-V^. (42) 



The explicit form is: 



24>-4> 2 + 2^Y / H i -Y t (2H i + Hfj - + 

i i \ i / 

8V 

+V{4>)- -==A s d -V<7(<£), (43) 

Q(j> 

where we have defined, by analogy with Eq. (6), 

a@) = -2%. (44) 

The new set of equations (9), (41), (43) is compatible with scale- factor 
duality for any V and a, as can be shown by rewriting the equations in terms 
of the shifted variables of Eq. (13). With such variables, Eqs. (9), (41), (43) 
become, respectively, 

f -J2 H ?- V = 2X t 1 ^P, (45) 

i 

H i -H l ^=\ d - 1 ^p i , (46) 
2% - f - £ Hf + V{4>) - ^ = A^ 1 e^ a®). (47) 

They are manifestly invariant under the generalized transformations 

(H^aT 1 , ~4>^~4>, P^P, Pi^-Pi, ct^ct, (48) 
preserving the shifted version of the dilaton-charge density a. 
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1.3 The pre-big bang scenario 

Let us come back to the cosmological applications of scale-factor duality. Even 
without using its non-local extensions, the duality symmetry of the equations 
allows introducing a "dual complement" of the standard cosmological solu- 
tions, and suggests new possible scenarios for the primordial evolution of our 
Universe. 

For a simple illustration of this possibility it will be enough to consider 
a homogeneous, isotropic and spatially flat metric background, sourced by a 
barotropic perfect fluid with equation of state vl P — l — const, with negligible 
dilaton charge. By imposing a = 0, V = 0, and assuming 7^0, one easily 
finds that Eqs. (45)-(47) are satisfied by the following particular exact solution 

°=UJ ' ~ p=poa ' ^ = -Tw ln UJ +const ' (49) 

where t > 0, and to, p are positive integration constants. In terms of the 
non-shifted variables: 



2-1 



to 



p = p a d(1+ ~<\ p = 7P, 



= 2 1 In ( — ) + const, t > 0. (50) 



l + d 7 2 V^o 

This solution, defined over the real positive semi-axis t > 0, describes a Uni- 
verse evolving from a past curvature singularity at t — > 0_ to an asymptoti- 
cally flat configuration at t — > +00. For 7 > we have a phase of decelerated 
expansion and decreasing curvature, 

a >0, a < 0, F<0, (51) 

typical of the standard cosmological scenario. Also, for a "realistic" equation 
of state with jd < 1, the dilaton turns out to be non- increasing (</> < 0); 
in particular, for a radiation fluid with 7 = 1/cZ, one recovers the radiation- 
dominated solution at constant dilaton, 

2 

P=- d , a= (^) 1+d ' P = Paa- {1+d \ cj> = const, (52) 

which is also an exact solution of the standard Einstein equations. 

Thanks to the symmetries of the string cosmology equations we can now 
obtain new, different solutions (which have no analogue in the context of the 
Einstein equations) by performing a time reflection t —* —t and, simultane- 
ously, a dual transformation defined by Eq. (18). Starting in particular from 
Eq. (50) we are lead to the background 
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• ^) , P = Paa d{1 7) , P=-7P, 

1 + d'Y / t \ 

0=-2 r^ ln + const, t < 0, (53) 

1 + efy 2 V t / 

which is still a particular exact solution of Eqs. (45)-(47). It is defined on the 
negative real semi-axis t < 0, and for 7 > it describes a phase of accelerated 
(i.e. inflationary) expansion and growing curvature: 

a > 0, a > 0, H > 0. (54) 

In this case the Universe evolves from an asimptotically flat initial configu- 
ration at t — > —00 towards a curvature singularity at t — > 0_. The dilaton 
is always growing (0 > 0) for t — > 0_, even if we consider the dual of the 
radiation-dominated solution (52): 

p = - d > a= {-T ) ' ^ = ^ oa < ^ = -rfTT ln l^ 

(55) 

This interesting property of the low-energy string-cosmology equations 
- i.e. the presence of an inflationary "partner" associated to any standard 
decelerated solution - is also valid in the absence of sources. Consider, for 
instance, Eqs. (45)-(47) with p = = p, and V = 0. In the isotropic limit we 
find the particular exact solution 



a = 



tV /V ~ d f r, \, ft 



to J V / V to 



= (Vd- l) In ( — J , t>0, (56) 



describing decelerated expansion and decreasing curvature. By applying the 
transformations (18) we are lead to the dual solution 



; \ -1/Vdi 



= - (Vd + l) In (-j-J , t<0, (57) 



to J ^ ' \ to 

describing accelerated expansion, growing curvature and growing dilaton. Ac- 
tually, both the vacuum and the fluid-dominated solutions can be obtained as 
asymptotic limits of the general exact solution of the system of equations (45)- 
(47) (for barotropic sources, with V = and a = 0), in the large-curvature 
and small-curvature limits, respectively [12, 13]. 

It is well known that the decelerated configurations, typical of the standard 
cosmological scenario, cannot be extended back in time without limits: the 
range of the time coordinate is bounded from below by the presence of the 
initial singularity (indeed, going back in time, the growth of H is unbounded, 
and the curvature blows up to infinity in a finite proper-time interval). 
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The standard scenario, however, is certainly incomplete because it ex- 
cludes inflation. The inclusion of inflation, on the other hand, modifies the 
behavior of the curvature scale: during a phase of "slow-roll" inflation [14], 
for instance, the background geometry can be approximately described by a 
de Sitter-like metric where H ~ 0, and in which the curvature tends to set- 
tle at a constant. One might think, therefore, that a complete (and realistic) 
cosmological scenario could avoid the initial singularity, replacing it with a 
primordial inflationary phase at constant curvature. 

Unfortunately, however, an epoch of accelerated expansion at constant 
curvature, described by the Einstein equations, and dominated by the poten- 
tial energy of some "inflaton" scalar field satisfying causality and weak-energy 
conditions, cannot be "past eternal", as proved in [15]. Thus, the conventional 
inflationary scenario mitigates the rapid growth of the curvature typical of the 
standard cosmological evolution, and shifts back in time the position of the 
initial singularity, without completely removing it, however (namely, without 
extending in a geodesically complete way the model, back in time, to infinity). 

If a constant-curvature phase is not appropriate to construct a regular 
model (fully extended over the whole temporal axis), the alternative we are 
left is a model in which the curvature, as we go back in time, after reaching a 
maximum, at some point, starts decreasing; in other words, a model in which 
the standard evolution is completed and complemented by a primordial phase 
with a specular behavior of H with respect to the standard one. Remarkably, 
this is exactly what can be obtained assuming that the cosmological evolution 
satisfies a principle of "self-duality" - i.e. assuming that the past evolution of 
our Universe is described by the "dual complement" of the present one [2, 22]. 

More precisely, if we consider a cosmological model satisfying (at least ap- 
proximately) the self-dual condition a(t) = a _1 (— f), such that the standard 
decelerated regime at t > smoothly evolves, back in time, into the acceler- 
ated partner at t < 0, we can then obtain a scenario in which the singularity 
is automatically regularized, and the initial evolution is automatically of the 
inflationary type. In such a context the big bang singularity is replaced by an 
epoch of high (but finite) curvature, characterizing the transition between the 
standard cosmological phase (H < 0) and its dual (H > 0): it comes natural, 
in such a context, to call "pre-big bang" the initial phase (f < 0) at growing 
curvature and growing dilaton, in contrast to the subsequent "post-big bang" 
phase (t > 0), describing the standard cosmological evolution. 

The dilaton, on the other hand, provides an exponential parametrization of 
the (tree-level) string coupling g s = exp(</>/2), controlling the relative strength 
of all (gravitational and gauge) interactions [3,4]. The principle of self-duality 
thus suggests that the Universe is lead to its present state after a long evolution 
started from an extremely simple almost trivial - configuration, character- 
ized by a nearly flat geometry and by a very small coupling parameter, 

H 2 ^0, cj> --oo, <£ = e*->0, (58) 
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the so-called string perturbative vacuum (see Fig. 2). In this case, the initial 
Universe is characterized by a regime of extremely low energies in which the 
"curvatures" (i.e. the field gradients) are small (X 2 H 2 <C 1, \ 2 <fi 2 < 1, ...), 
the couplings are weak (g 2 <§; 1), and the background dynamics can be appro- 
priately described by the lowest-order string effective action, at tree-level in 
the a' and quantum loop expansion (also in agreement with the hypothesis of 
"asymptotic past triviality" [16]). We can talk of "birth of the Universe from 
the string perturbative vacuum" , as also pointed out in a quantum cosmology 
context (see e.g [17, 18]. 
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Fig. 2. Qualitative time-evolution of the curvature scale (upper panel) and of the 
string coupling (lower panel), for a typical self-dual background which smoothly 
interpolates between the pre-big bang and the post-big bang phase, starting from 
the string perturbative vacuum. 



This picture is in remarkable contrast with the standard (even inflationary) 
picture in which the Universe starts evolving from a highly-curved geomet- 
ric state: the more we go back in time, in that context, the more we enter 
a Planckian and (possibly) trans-Planckian [19] non-perturbative regime of 
ultra-high energies, requiring the full inclusion of quantum gravity effects, to 
all orders, for a correct description. 
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The principle of self-duality, on the contrary, suggests a picture in which 
the more we go back in time (after crossing the epoch of maximal curvature) , 
the more we approach a flat, cold and vacuum configuration (strongly reminis- 
cent of the "biblical" scenario quoted at the beginning of Sect. 1), which can be 
appropriately described by the classical background equations obtained from 
the action (1). Quantum effects, in the form of higher-curvature and higher- 
loop contributions, are expected to become important only towards the end of 
the pre-big bang phase, when the background approaches the string scale at 
t — > 0_. Actually, all studies performed so far have shown that such correc- 
tions must become dominant, eventually, in order to stop the growth of the 
curvature [20] and possibly trigger a smooth transition to the post-big bang 
regime [21]. 

1.4 A smooth "bounce" 

The lowest-order string effective action can appropriately describe the phase 
of primordial background evolution typical of the pre-big bang scenario, but 
not the transition to the standard decelerated regime occurring at high cur- 
vatures and strong coupling, and requiring the introduction of higher-order 
corrections. Referring the reader to the existing literature for a detailed re- 
view of the transition models studied so far (see for instance [22]), we shall 
present here only two simple phenomenological examples, by applying, to this 
purpose, the formalism introduced in subsection 1.2 (and Appendix A). In 
these examples, in fact, the bouncing transition is induced by the presence of 
a non-local effective potential V{<f>), expected to simulate the backreaction of 
the quantum loop corrections in higher-dimensional manifolds with compact 
spatial sections [11]. 

The first example is based on a potential which, in the homogeneous limit, 
takes the form 



and which may thus perturbatively interpreted as a four-loop potential. With 
this potential, the duality-invariant equations (45)-(47), in vacuum (p = p = 
<T = 0), and in the isotropic limit, are solved by the particular exact solution 



v{<j>) 




Vb > 



(59) 



[18]: 




<fi = In 



+ const, 



(60) 
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where to and a>o are positive integration constants. This regular "bouncing" 
solution is exactly self-dual - as it satisfies a(t)/a = a /a(—t) - and is char- 
acterized by a bounded, "bell-like" shape of the curvature and of the dila- 
ton kinetic energy (see Fig. 3). The solution smoothly interpolates between 
the pre-and post-big bang vacuum solutions (57), (56) (corresponding to the 
dashed curves of Fig. 3), which are recovered in the asymptotic limits t — > — oo 
and t — > -(-oo, respectively. The bounce of the curvature, and the smooth tran- 
sition between the two branches of the low-energy solutions, is induced and 
controlled by the potential (59) which dominates the background evolution 
in the high-curvature limit \t\ — > 0, and which becomes rapidly negligible as 
t — » ±oo, as illustrated in Fig. 3. 



Fig. 3. Plot of the curvature, of the dilaton kinetic energy, and of the potential 
V((f>), for the bouncing solution (60). The dashed curves represent the (singular) 
vacuum solutions (56), (57), obtained with V = 0. All curves are plotted for to = 1, 
Vo = 1, and d — 3. 

It should be noted that in this solution the dilaton keeps growing, mono- 
tonically, even in the limit t — > +oo. In more realistic examples, however, 
such a growth is expected to be damped by the interaction with the mat- 
ter/radiation post-big bang sources [23], and/or by the action of a suitable 
non-perturbative potential appearing in the strong coupling regime. 

The second example of bounce is based on a general integration of the 
duality-invariant equations (45)— (47), in the presence of isotropic fluid sources 
with a — and of a two-loop (non-local) potential which in the homogeneous 
limit takes the form 




t 



v{4>) 




V >0. 



(61) 



In this case the equations can be integrated exactly not only for barotropic 
equations of state (p/p = 7 = const), but also for any ratio p/p which is an 
integrable function of an appropriately defined time-like parameter [2] . 
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An interesting example (motivated by the study of the equation of state 
of a string gas in rolling backgrounds [24]) is the case in which p/ p smoothly 
evolves from the value 7 = — 1/d at t = —00 to the value 7 = 1/d at t = +00, 
thus connecting the radiation equation of state to its dual partner, according 
to the law: 



P 1 



(62) 



Here x\ is an arbitrary integration constant, and a; is a (dimensionless) time- 

dx L 



like coordinate defined by 



at = 2 P ' 



(63) 



where L is a constant with dimensions of length (we are using units in which 
2Af _1 = 1, so that [p] = L~ 2 ). Using Eqs. (61)-(63), and choosing a simpli- 
fying set of integration constants (appropriate to the pedagogical purpose of 
this paper), we can then obtain the following particular exact solution [2], 

2/{d-l) 



a = do x 



e^ = a n e T 



2d/(d-l) 



1 + 



pe* 
pe 4 



d- 1 
d-l 



^x 2 + x\ J 

e 2 ^ (x 2 + x\)- {d+1)/{d - l \ 

200^^2 , „2^-(3d+l)/2(d-l) 



( ^2^2" e " V " x \ x + X l 



(64) 



where ao and (pg are integration constants. The smooth and bouncing behavior 
of this solution is illustrated in Fig. 4. 

The above solution is self-dual, in the sense that <j){x) = <f>{—x), p(x) = 
~p(—x), and 



(65) 



a(x) 




a(—x) 


_aox 1 




n ^V(d-l) 

_aox 1 



(with an appropriate choice of the integration constant a it is always possible 
to set to 1 the fixed point of scale- factor inversion). The solution satisfies, 
asymptotically, 



-00 



+00 



r 2/(d-l) 



dx 

at' 

dx 



(66) 



Re-expressing a, <fr, p, p, in the asymptotic limits x — > ±00 in terms of the 
cosmic time t, we can check that this solution smoothly interpolates between 
the pre-big bang configuration (55) describing accelerated expansion, grow- 
ing dilaton, negative pressure, and the final post-big bang configuration (52), 
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Fig. 4. Plot of the curvature, of the string coupling, of the effective energy density 
and of the effective pressure for the self-dual solution (64). The curves are plotted 
for d = 3, L = 1, Xi = 1, 4>o = 0, and ao = exp(— 2/3). 



describing the radiation-dominated state with frozen dilaton and decelerated 
expansion. As in the previous case the smoothing out of the tree-level sin- 
gularity, and the appearance of bouncing transition, is a consequence of the 
effective potential (61). 



1.5 Cosmological perturbations 

The phase of pre-big bang evolution, being accelerated, can amplify the quan- 
tum fluctuations of the metric tensor (and of other background fields) just like 
any other type of inflationary evolution. However, because of the kinematic 
properties of pre-big bang inflation (associated to the shrinking of the Hubble 
horizon -ff -1 ), the spectral distribution of the metric fluctuations, after their 
amplification, tends to grow with frequency [25]. This peculiar aspect of the 
spectrum may be regarded as representing both an advantage and a difficulty 
of pre-big bang models with respect to other models of inflation. 

The advantage is of phcnomenological nature, and refers to the transverse 
and traceless tensor part of the metric fluctuations. Their amplification leads 
to the formation of a stochastic background of relic gravitational waves whose 
spectral energy density, fi g , grows wih frequency: 



(4i 



6>0, 



to < u±. 



(67) 



Here Mp = 8itG = Ap 1 is the Planck mass, Hi ~ M s — X^ 1 the inflation- 
radiation transition scale (expected to be controlled by the string mass scale 
Ms), Q r — Pr/pc is the fraction of critical energy density in radiation, loi is 
the ultraviolet cut-off (i.e. the maximal amplified frequency) of the spectrum, 
and 6 a model-dependent parameter depending on the background kinematics 
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[25, 26, 27] (see also the contribution of Buonanno and Ungarelli [28] to this 
volume) . 

Thanks to the growth of the spectrum, the cosmic graviton background 
present today as a relic of the inflationary epoch is higher at higher frequencies 
(in particular, higher than the backgrounds predicted by conventional models 
of inflation) , and thus more easily detectable by current gravitational antennas 
(see e.g. [22]). Conversely, however, the spectrum is strongly suppressed in 
the low-frequency regime: we should thus expect, in particular, a negligible 
contribution of tensor metric perturbations to the observed CMB anisotropy 
on large scales (as in the case of the ekpyrotic [29] and "new ekpyrotic [30] 
scenarios where, however, the gravitational background is expected to be low 
even in the low-frequency regime [31]). It may be stressed, in this connection, 
that the possible absence of tensor contributions at large scales emerging 
from (planned) future measurements of the CMB polarization (such as those 
of WMAP, PLANCK), in combination with a positive signal possibly detected 
at high frequency by the next generation of gravitational antennas (such as 
LIGO/VIRGO, LISA, BBO, DECIGO), could represent a strong experimental 
signal in favor of models of pre-big bang inflation (see e.g. [32]). 

The difficulties associated to a growing spectrum refer to the scalar part of 
the metric perturbations. In fact, a growing scalar spectrum cannot account 
for the observed peak structure of the temperature anisotropics of the CMB 
radiation, which requires, instead, a nearly flat (or "scale-invariant" ) primor- 
dial distribution: Q a {w) <~ w" s_1 , with n s w 1. There are two possible ways 
out of this problem. 

A first possibility relies on the growth of the dilaton - and thus of the string 
coupling gl = exp <fi - during the phase of pre-big bang inflation. Even starting 
at weak coupling, a pre-big bang background unavoidably evolves towards the 
strong coupling regime gr s ~ 1. If the bounce is not immediate then the Uni- 
verse, before the transition to the standard regime, enters a strong-coupling 
phase where higher-dimensional extended objects like Dirichlet branes and 
antibranes [4] (whose tension is proportional to the inverse of the string cou- 
pling) become light, and can be copiously produced [33]. The cosmic evolution 
may become dominated by the presence of these higher-dimensional sources 
[34] and, in that context, a phase of conventional slow-roll inflation can be 
triggered by the interaction (and eventual collision) of a brane-antibrane pair 
[35] (see also the contribution of Tye [36] to this volume) . This new inflation- 
ary regime may efficiently dilute all pre-existing inhomogeneities and generate 
a new spectrum of scale-invariant, adiabatic scalar perturbations, as required 
for a successful explanation of the observed anisotropy. This may resolve the 
incompatibility between a (growing) spectrum of pre-big bang perturbations 
and present large-scale observations. 

There is, however, a second possibility which avoids introducing additional 
inflationary epochs besides the initial dilaton-dominated one, and which is 
based on the so-called "curvaton mechanism" [37]. According to this mecha- 
nism the (flat, adiabatic) spectrum of scalar metric perturbations, responsible 
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for the observed anisotropies, is not produced during the primordial evolution: 
instead, it is the outcome of the post-inflationary decay of a massive scalar 
field (the curvaton), whose quantum fluctuations are amplified during inflation 
with a nearly flat spectrum, and are converted into curvature perturbations 
after its decay. In the context of the pre-big bang scenario the role of the 
curvaton is possibly played by the Kalb-Ramond axion a [38], associated - 
by space-time duality - to the four components of the NS-NS two-form B^ v 
present in the massless multiplet of the string spectrum. 

For a brief discussion of this possibility we should explain, first if all, 
why axion fluctuations can be amplified by pre-big bang inflation with a flat 
spectrum [39], unlike metric fluctuations. The reason is that the slope of the 
spectrum is directly related to kinematic behavior of the effective "pump field" 
responsible for the amplification, and that metric and axion fluctuations have 
different pump fields, even in the same given background. 

In order to clarify this point let us complete the low-energy action (1) by 
adding the contribution of the antisymmetric field B^ v , considering (for sim- 
plicity) a model already dimensionally reduced to four space-time dimensions: 

= d^B va + d v B ail _ + d a B^ u . (68) 

In the absence of sources the equations of motion for B^ are automatically 
satisfied by introducing the "dual" axion field a, such that 

H^ a = -j= e^dpa, (69) 

and the last term of the action (68) can be replaced by 

S=^J ^x^-ge^Vaf. (70) 

Perturbing the metric and the axion field, 

9y. v — > g^v + h^, a —> a + Sa, (71) 

around a homogeneous, conformally flat metric background, using the con- 
formal time coordinate r] (such that dt = adt]), and applying the standard 
formalism of linear cosmological perturbations (see e.g. [40]), we obtain for 
tensor metric and axion fluctuations, respectively, the following quadratic ac- 
tions: 

Sh = \J d 3 xdrjz 2 h (rj)(h' 2 + hV 2 h), 

*. - a p-^/ 2 (t?:\ 
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S* = \j d 3 xd V z 2 J V ) (5a' 2 + 5aV 2 5a) 



Za = 



V2A S 



>/2 



(73) 



Here h is one of the two physical polarization states of tensor perturbations, 
the primes denote differentiation with respect to 77, and V 2 is the flat-space 
Laplace operator, V 2 = S^didj. The variation of these actions with respect 
to h and 5a leads to the equations of motion, which can be written in terms 
of the canonical variables u = (hzh) and v — (Saz^) as follows: 



(hz h y - 

(5az a )" 



=*■ ) (hz h ) 

Zh, 



0, 



V 2 + ^ 



(5az a ) = 0. 



(74) 
(75) 



The canonical equations are the same for u and v, but the pump fields, Zh 
and z ai are different. 

Consider, for instance, the axion equation (75), and recall that during 
inflation the accelerated evolution of the pump field can be parametrized as 
a power-law evolution in the negative range of the conformal-time parameter 
[22, 40], i.e. 



OO < T] < 0, 



(76) 



where rji > is some appropriate reference time-scale. Expanding in Fourier 
modes, Eq. (75) becomes a Bessel equation for the mode Vk, 



4 + 



k 2 



(»l 1/4) 



vk = 0, 



1 



(77) 



and its general solution can be conveniently written as a combination of first- 
kind and second- kind Hankcl functions [41], of argument krj and index v a , as 
follows: 



v k = (-??) !/ 2 A + {k)H^{krf) + A_{k)H^{krf) 



(78) 



We shall now canonically normalize this general solution by imposing that 
the initial state of the fluctuations corresponds to a spectrum of quantum 
vacuum fluctuations [22, 40]. More explicitly, we shall require that the mode 
Vk, on the initial spatial hypersurface at r\ — > — 00, may represent freely oscil- 
lating, positive frequency modes satisfying the canonical normalization 



from which 



V k V k - V k V k = 2, 



-ikr] 



v k 



2k 



-00 



(79) 



(80) 
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(modulo an arbitrary phase). Using the large argument limit of the Hankel 
functions [41], 



HP ikr,) = ^ e"— , (k v ) = ^ (81) 

(e„ = — 7r/4 — vk/2), we obtain A + = ^/7r/4 and A- = 0. The normalized 
exact solution for the the axion fluctuations 5<Tk can be finally written as 

where 6^ is an arbitrary phase determined by the choice of the initial condi- 
tions. 

In order to determine the spectrum of the fluctuations after their inflation- 
ary amplification we must then consider the limit — > 0_ , in which \kt]\ <C 1 
and the amplitude of the mode k is stretched "outside the horizon" . We can 
use, to this purpose, the small argument limit of the Hankel functions [41], 
which reads (for v ^ 0), 

HPikv) = V l{kriT - iq v {kr,)- v + ... (83) 

where q v and p v are complex (^-dependent) coefficients (for v = there are 
additional logarithmic corrections). We obtain, in this limit, 



° k ~ z a ~* Mp \ 2 J 



»?1 



(84) 



The cases we are interested here are limited to "conventional" inflationary 
backgrounds with a a < 1/2, i.e. v a > (sec [32] for a detailed discussion of 
all possibilities). For such backgrounds the time-dependence of 8a k tends to 
disappear as 77 — > 0_, the fluctuations become frozen, asymptotically, and their 
(dimensionless) spectral amplitude k 3 \5<Jk\ 2 , controlling the typical amplitude 
of the perturbations on a comoving length scale r — k~ 1 [40], has the following 
fc-dependencc: 

k 3 \6a k \ 2 ~ k 3 - 2l/ ° = k 2+2a °. (85) 

This result also holds in the limiting case a a — 1/2 with the only addition of 
a mild logarithmic correction [26, 27], i.e. fc 3 |<5cr/c| 2 ~ fc 3 ln 2 (fc?7i). 

The above calculations can be exactly repeated, in the same form, for the 
tensor perturbation variable, starting from Eq. (74): the resulting spectrum 
is formally the same, 

k 3 \h k \ 2 ~ k 3 - 2 » h =k 2+2a \ (86) 

with the difference that the spectral slope is now determined by the power 
ah, controlling the evolution of the tensor pump field Zh through an equation 
analogous to Eq. (76). 
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We are now in the position of discussing the possible pre-big bang pro- 
duction of a fiat spectrum of axion fluctuations, even if the associated metric 
fluctuations are amplified (in the same background) with a growing spectrum. 
Let us consider, to this purpose, an exact anisotropic solution of the string 
cosmology equations (9)— (11), in vacuum, and without dilaton potential. The 
solution describes a phase of pre-big bang inflation characterized by the ac- 
celerated (isotropic) expansion of three spatial dimensions, with scale factor 
a(rj), and by the accelerated contraction of n "internal" spatial dimensions, 
with scale factors 6^(77), i = 1, . . . ,n. In conformal time, such a solution can 
be parametrized for 77 — > 0_ as [13, 22] 

n /3o/(l-/3o) / \/3i/(l-A>) 



m/ Km 

^ = E.ft+y.-i ,/,y (87) 



i-A) V m. 

where the constant coefficients (3 , Pi satisfy the Kasner-like condition 

> 2 + 3/3 2 = l, (88) 



E^ 2 



and 4>4, +n is the higher-dimensional dilaton appearing in the full (4 + n)- 
dimensional effective action. The four-dimensional dilaton cf) is related to 4>i+n 

by 

e~* = V n e"^+" = e-* 4 +» [] h, (89) 



namely by 



= ^ + „-^ln6 i = ^-ilnf-^-). (90) 



Let us compute, for this background, the kinematic powers and a a 
controlling the evolution of the pump fields (72), (73): 

a ^w^h- (92) 

It follows, according to Eq. (86), that the spectrum of tensor (as well as of 
scalar) metric perturbations is always characterized by a slope which is cubic 
(modulo log corrections) [13, 26, 27], and which is also "universal", in the 
sense that it is insensitive to the background parameters (n, /?o, Pi)- For the 
axion fluctuations, on the contrary, we find from Eq. (85) that the spectral 
slope is strongly dependent on such parameters, and that a scale-invariant 
spectrum with 2 + 2a a = is allowed, in particular, provided /3 = —1/3. 

We may note, in the special case in which the background is fully isotropic 
and expanding (i.e., /3q = Pi < 0), that the Kasner condition (88) implies 
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(3q = —l/Vd, so that a scale-invariant spectrum corresponds to d = 9, i.e. 
just to the number of spatial dimensions determined by critical superstring 
theory [3, 4]. 

In the less special case in which the spatial geometry can be factorized as 
the product of a 3-dimcnsional and a n-dimcnsional isotropic subspaces we 
have, instead, f3i — (3 ^ /3 , and 3j3 2 + n[3 2 — 1. The spectral slope, in this 
case, can be expressed in terms of the parameter 

-i(£)(sr-& 

controlling the relative time-evolution of the proper volumes of the internal 
and external spaces. Eliminating [3 in terms of /3 through the Kasner con- 
dition, and replacing (3q with r in Eq. (92), one can then parametrize the 
deviations from a flat axion spectrum as the relative shrinking or expansion 
of the two subspaces [42] . 

Given a sufficiently flat spectrum of axion fluctuations, amplified by the 
phase of pre-big bang inflation, we are then lead to a post-big bang configura- 
tion which is initially characterized (at some given time scale %) by a primor- 
dial sea of "isocurvature" scalar perturbations, dominated on super- horizon 
scales by the axion fluctuations 6a (the metric fluctuations are subdominant 
on such large scales, being strongly suppressed by the steep slope of their 
spectrum). The axion can play the role of the curvaton provided that the ini- 
tial configuration, besides containing the initial fluctuations 6a i, also contains 
a non- vanishing axion background, <r, ^ 0, whose energy density p a - even if 
subdominant - is initially determined by an appropriate potential (possibly 
approximated by V a ~ m 2 a 2 ). In that case the background evolution, after 
an initial slow-roll regime, leads to a phase where the axion background starts 
oscillating with proper frequency m, at a curvature scale H ~ m, simulating 
a dust fluid (p a ~ a~ 3 ) which may become dominant with respect to the 
radiation fluid, and eventually decay at the typical scale H ~ Apm 3 . 

In such a type of background the axions fluctuations 6a become linearly 
coupled to scalar metric perturbations, and may act as sources for the so- 
called Bardeen potential <P. New metric perturbations can then be generated, 
starting from ^(rji) — 0, with the same spectral slope as the axion one, and 
with a spectral amplitude not smaller, in general, than the axion amplitude. 
Referring to the literature for a detailed computation [37, 38], we shall recall 
here that the final spectrum (after the axion decay) of the super-horizon 
Bardeen potential is related to the initial axion perturbations by 

\&k\ = A P /(a;) \6a k (-m)\ , 

f{(7i) = ci-^ 2 - + c 2 + c 3 X P ai (94) 

(the Ap factors are due to the canonical normalization of the axion field and 
of its fluctuations). Here is the initial amplitude of the axion background, 
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fl a <~ 1 is the axion fraction of critical density at the axion decay epoch, and 
c ij c 2j C3 are dimensionlcss numbers of order one {fi a cannot be much smaller 
than one, to avoid a too strong "non-Gaussianity" of the spectrum" [43]). 
Thanks to its structure, the "form factor" f{<Ji) has a minimum of order one 
around Apc^ ~ 1. A (nearly) scale- invariant axion spectrum thus reproduces 
a (nearly) scale- invariant spectrum of scalar metric perturbations. 

As discussed in the literature, a curvaton-induced spectrum of scalar met- 
ric perturbations provides the right "adiabatic" initial conditions for repro- 
ducing the observed temperature anisotropies of the CMB radiation, exactly 
as in the case of the slow-roll scenario. The only difference is the "indirect" 
(i.e., post-inflationary) production of the scalar spectrum, triggered by the 
presence of a non-vanishing axion background. It must be stressed, however, 
that the direct connection (94) with the axion spectrum of primordial origin 
gives us the possibility of extracting, from present CMB observations, impor- 
tant constraints on the parameters of pre-big bang models of inflation [38] . 

In particular, using the experimental normalization of the anisotropy spec- 
trum, and the direct relation between the pre-big bang inflation scale H\ and 
the string scale M s , one can speculate about the possibility of "weighing the 
string mass with the CMB data" [44] . Another application concerns the slope 
of the scalar perturbation spectrum which, according to most recent WMAP 
results [45], is given by 

n s = 3 + 2a = 0.951±g$g. (95) 
Using Eq. (92), and the Kasner condition (88), one obtains 

(3o- -0.355, ^/? 4 2 ~0.62. (96) 

i 

With d = 9 dynamical dimensions this result seems to point out the existence 
of a small anisotropy between the kinematics of the external and internal 
spaces during pre-big bang inflation (a fully isotropic expansion would corre- 
spond, in fact, to (3 = -I/V9 ~ -0.33 and = 6/9 ~ 0.66). It should 
be noted, however, that other interpretations of the data are also possible. 
For instance, the result (95) is also compatible with /? = — 1/V8 ~ —0.3535, 
describing the isotropic expansion of d = 8 spatial dimensions! Incidentally 
the number (and the kinematics) of the extra spatial dimensions play a crucial 
role also in the possible production of primordial "seeds" for the large-scale 
magnetic fields [46]. 

It should be mentioned, finally, a possible non-Gaussian "contamination" 
of the statistical properties of the anisotropy spectrum, possibly present in 
curvaton models with fi a <C 1 [43] (see Eq. (94)). A possible detection of 
non-Gaussianity, in future CMB measurements, could provide support to the 
curvaton mechanism, and could be used for a direct discrimination between 
this scenario and other, more standard scenarios based on slow-roll inflation. 
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2 The relic dilaton background 

The accelerated evolution of the Universe, during the phase of pre-big bang in- 
flation, amplifies the quantum fluctuations of all fields present in the string ef- 
fective action: thus, in particular, it amplifies the dilaton fluctuations, 8<j> = X- 
The formation of a stochastic background of relic gravitational waves, asso- 
ciated to the amplification of the tensor part of metric fluctuations, is thus 
accompanied by the simultaneous formation of a comic background of relic 
dilatons [47], whose primordial (high-energy) spectral distribution tends to 
follow that of tensor metric perturbations [13]. 

There is, however, a possible important difference in the present intensity 
of the two cosmic backgrounds, due to the fact that dilatons - unlike gravi- 
tons - could become massive in the course of the standard (post-inflationary) 
evolution. Actually, dilatons must become massive if they are non-universally 
coupled to ordinary matter with gravitational strength (or higher) [48, 49], 
to avoid the presence of long-range scalar forces which are excluded by the 
standard gravitational phenomenology (in particular, by the high-precision 
tests of the equivalence principle). The induced mass may drastically modify 
the amplitude and the slope of the dilaton spectrum, in the frequency band 
associated to its non-relativistic sector. 

For a simple illustration of the effects of the mass on the spectrum we 
will consider here the model of vacuum, dilaton-dominated pre-big bang back- 
ground described by Eq. (57), smoothly joined at 77 = —771 < to the standard 
radiation-dominated background with frozen dilaton, described by Eq. (52) 
(we shall work in d — 3 spatial dimensions) . Perturbing the background equa- 
tions [13] one finds, in this case, that the dilaton pump field is the same field 
Zh ~ aexp(— <f>/2) governing the amplification of metric fluctuations. Taking 
into account a possible mass contribution, m 2 = d 2 V/d(f> 2 , one then obtains 
for the Fourier modes Xfe the canonical equation: 



During the initial pre-big bang regime the potential is negligible (m 2 = 0), 
and the canonically normalized solution for Xfe is that of Eq. (82) (with v a 
replaced by Vh). In the subsequent radiation-dominated era <p stabilizes to a 
constant, so that Zh ~ a ~ r\ and the effective potential z'^/zh is vanishing. 
Assuming that the dilaton mass is small enough in string units, and consid- 
ering the high-frequency sector of the spectrum, associated to the relativistic 
modes of proper momentum p = (k/a) 3> m, we can neglect also the mass 
term of Eq. (97), to obtain the general solution 




(97) 



1 



[c + (fc)c- ^fer '+c_(fc)e ^ ' ^, '] 



(98) 



Xk = 




1] > -T}1- 



Matching \ an d x' with the pre-big bang solution (82) at rji, for super-horizon 
modes with (fcryi) Cl, we are lead to 
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c±(fc) = ±c(k)e* ik « 1 , \c(k)\ ~ {krjx)-^- 1 ' 2 (99) 

(modulo numerical factors with modulus of order one). Thus, at large times 
V > Vi, 

c(k) 

Xfc ~ sinfcr?. (100) 
a\Jk 

The spectral energy density for the relativistic sector of the dilaton back- 
ground, in the radiation era, in then determined by 

5)W~(SK±)~-/®~ ™ 

where k\ <~ 77 is the high-frequency cut-off scale. In units of critical energy 
density, p c = 3M|i7 2 , 

where we have defined the (model-dependent) slope parameter 6 = 3—2^ > 0, 
and we have introduced the (time-dependent) proper momentum associated 
to the cut-off scale, pi = k\/a — H\a\/a, determined by the background 
curvature scale H\ at the end of inflation. In general, (Hi/H) 2 (ai/a) 4: = 
p r (t)/p c (t) = f2 r (t), and we may thus conclude that the relativistic sector of 
the dilaton spectrum, in the radiation era, is exactly the same as the spec- 
trum of tensor metric perturbations (see Eq. (67)), in the same model of 
background. 

However, even if the mass is small, and initially negligible, the proper mo- 
mentum p = k/a{t) is continuously red-shifted with respect to to during the 
subsequent cosmological evolution, so that all modes tend to become non- 
relativistic, p < to. For non-relativistic modes the solution (98) is no longer 
valid, and the correct spectrum must refer to the exact solutions of Eq. (97) 
with to ^ 0. In the radiation era such a solution can be given in terms of the 
Weber cylinder functions [50], and one finds that the non-relativistic sector 
of the spectrum splits into two branches, with different slopes: a first branch 
of modes becoming non-relativistic at a time scale t nr when they are already 
inside the horizon, with proper momentum p such that p(t nr ) <~ m ^> H(t nr ); 
and a second branch of modes becoming non-relativistic when they are still 
outside the horizon, with p(t nr ) ~ m <C H(t nr ). The two branches are sepa- 
rated by the momentum scale p m of the mode becoming non-relativistic just 
at the time of horizon crossing, i.e. p(t nr ) = to = H(t nr ), and thus related to 
the cut-of scale p\ by 



/ tt \ 1/2 / \ 1/2 

p m TOa nr m Hi \ I jn 



pi Hiai Hi \H m J \Hi 



(103) 
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Without applying to the explicit form of the massive solutions of Eq. (97) , 
a quick estimate of the non-relativistic spectrum can be obtained [51] by 
noting that, if p nr > H(t nr ), the number of produced dilatons is the same as 
in the relativistic case, and the only effect of the non-relativistic transition is 
a rescaling of the energy density, i.e. 

flf - ^ = (j) ^\ (104) 

For this branch of the spectrum we then obtain, from Eq. (102), 




p m <p<m. (105) 



In the case p nr < H nr , on the contrary, the slope of the spectrum - 
determined by the background kinematics at the time of horizon exit - has to 
be the same as that of the relativistic sector, while the time-dependence has 
to be the non-relativistic one (p x ~ a~ 3 ) of Eq. (105). Continuity with the 
branch (105) at p = p m then gives 




Pcq < P < Pm- (106) 



The lower limit p e q < P has been inserted here to recall that we are ne- 
glecting the effects of the transition to the matter-dominate phase, i.e. we 
are considering modes re-entering the horizon during the radiation era, with 
p > p cq — H cq <~ 10~ 27 eV. We should recall, also, that the spectrum has been 
computed in a radiation-dominated background, and thus is valid, strictly 
speaking, only for t > i oq . 

The three branches (102), (105), (106) describe the spectrum (between 
p cq and pi) of primordial dilatons produced in the simple example of "mini- 
mal" pre-big bang model that we have considered. We refer to the literature 
for a more detailed computation, for a discussion of its transmission to the 
present epoch to, and for the possible modifications induced by generalized 
background evolutions (see e.g. [32]). For the pedagogical purpose of this pa- 
per this example provides a sufficiently clear illustration of the effects of the 
mass on the spectrum: in particular, it clearly illustrates the enhancement 
produced at lower frequencies because of the reduced spectral slope of the 
branch (105), which may become even decreasing if 5 < 1 (see Fig. 5). 

In such a context one is naturally lead to investigate whether this enhanced 
intensity might favor the detection of a non-relativistic dilaton background, 
with respect to other, relativistic types of cosmic radiation (such as the relic 
graviton background). 

2.1 Light but non-relativistic dilatons 

For a phenomcnological discussion of this possibility we must start with two 
important assumptions. The first is that the produced dilaton are light enough 
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to have survived until the present epoch. Supposing that massive dilatons have 
dominant decay mode into radiation (e.g., two photons), with gravitational 
coupling strength, i.e. with a decay rate r ~ A p to 3 , it follows that the primor- 
dial graviton background is still "alive" in the present Universe (characterized 
by the time-scale Hq 1 ) provided Hq 1 < r~ x , i.e. 



The second assumption we need is that the total energy density of the 
dilaton background, integrated over all modes, turns out to be dominated by 
its non-relativistic sector. Only in this case we can evade the stringent bound 
imposed by the nucleosynthesis, which applies to the relativistic part of any 
cosmic background of primordial origin. 

The energy density of a relativistic background, in fact, evolves in time like 
the radiation energy density, p rel / p ra d — f2 rel / f2 ra d = const: the present value 
of their ratio is thus the same as the value of the ratio at the nucleosynthesis 
epoch. To avoid disturbing the nuclear processes occurring at that epoch, 
on the other hand, one must require that Q rel / Q ra d < 0.1 [52]. Using the 
present value of fl ra d one is lead then to the constraint (l rel (to) ^ 5 x 10 -6 , 
which imposes a severe constraint on all relativistic primordial backgrounds. 
In particular, it imposes an upper limit on the peak value of the graviton 
background produced in models of pre-big bang inflation, thus determining 
the minimal level of sensitivity required for its detection [22]. 

The energy density of a non-relativistic background, on the contrary, 
evolves like the dark-matter density, and grows in time with respect to the 
radiation background: H nr / ' Q r ad ~ a. As a consequence, the value of i?" r can 
be very large today, even if negligible at the nucleosynthesis epoch. The only 
constraint we must apply, in this case, is the critical density bound, 



to be imposed at any time t, to avoid a Universe over-dominated by such a 
cosmic background of dust matter. Here h ~ 0.73 is the present value of the 
Hubble parameter Hq in units of 100 km s _1 Mpc. 

For the dilaton spectrum of Eqs. (102)-(106) there are, in particular, two 
different cases in which the total energy density is dominated by the non- 
relativistic modes. A first (obvious) possibility is the case in which all modes 
of the spectrum are presently non-relativistic, namely p\ (to) < to (in this case 
the branch (105) extends from p m to p\). This implies, however, that 



to < 10 2 MeV. 



(107) 




(108) 



Hid! 




1/2 



) 



2/3 



TO > 



(109) 
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For a typical string- inflation scale, Hi ~ M a , we obtain a lower limit on 
m which is well compatible with the upper limit (107), but which requires 
mass values too high to be compatible with the sensitivity band of present 
gravitational antennas (see Subsection 2.2). 

The second (more interesting) possibility is the case in which m < pi(to), 
but the parameter 5 is smaller than one, and the slope is flat enough, so that 
the spectrum is peaked not at pi but at p m = p\{m/ 'Hi) 1 ' 2 (see Fig. 5). 
In that case the momentum integral (108) is dominated by the peak value 
f2 x (p m ), and the critical density bound can be approximated by the condition 
Sl x (p m ,to) < 1. Using Eq. (105), and noting that in the matter-dominated 
era (t > t cq ) the value of the non-relativistic spectrum keeps frozen at the 
equality value f2 x (t eq ), we are lead to the condition fl x (t cq ,p m ) < 1, which 
implies 

m< (H cq A4Ht 4 ) 1/(5+1) . (HO) 
For Hi ~ M s , and (5 — * 0, this bound can be saturated by masses as small as 

m~ff eq (^|y~10- 23 eV. (Ill) 
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Fig. 5. Example of dilaton spectrum dominated by the non-relativistic sector. The 
spectrum is peaked at p = p m , and the slope parameter satisfies the condition S < 1. 

It is quite possible, therefore, to have a dilaton mass small enough to 
fall within the sensitivity range of present gravitational detectors, even if the 
energy density of the dilaton background is dominated by non-relativistic 
modes (thus evading the relativistic upper bound fl rel < 10 -6 ), and even 
if the background intensity is large enough to saturate the critical density 
bound, Q x ~ 1. 

So small mass values, however, are necessarily associated with long-range 
dilaton forces: in particular, if the mass satisfies the condition m < pi(to) ~ 
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(M S /M P ) 1 / 2 10 -4 eV (as in the example illustrated in Fig. 5), the correspond- 
ing force has a range exceeding the centimeter. This might imply macroscopic 
violations of the equivalence principle (due to the non-universality of the dila- 
ton coupling [48]), and macroscopic deviations from the standard Newtonian 
form of the low-energy gravitational interactions (which seem to be excluded, 
however, by present experimental results [53, 54]). 

We should recall, in fact, that in the presence of long-range dilaton fields 
the motion of a macroscopic test body with nonzero dilaton charge is no 
longer described by a geodesies. There are forces on the test body due to 
the gradients of the dilaton field, according to the generalized conservation 
equation 

V„T M v = |v„0, (112) 

following from the application of the contracted Bianchi identity to the gravi- 
dilaton equations (3) and (5). The integration of this conservation equation 
over a (space-like) t = const hypersurface then gives, in the point-particle (or 
monopole) approximation, the non-geodesic equation of motion [55] 

— + r a ^u a u? = ,V^, (113) 

where q is a dimensionless ratio representing the relative intensity of scalar to 
tensor forces (i.e., the effective dilaton charge per unit of gravitational mass 
of the test body) . 

For the fundamental components of macroscopic matter, such as quark and 
lepton fields, the value of q (or of the charge density a) is to be determined 
from an effective action which includes all relevant dilaton loop corrections 
[48, 13], and which is of the form 



z r (<i>)r- z+mvt? -v(<p) 



+ ZimV^) 2 - MfzU^f 



(114) 



Here we have used, for simplicity, a scalar model of matter fields tpi, and we 
have called Z the dilaton "form factors" arising from the loop corrections. 
The effective dilaton charge, therefore, turns out to be frame-dependent (the 
charge q appearing in Eq. (113), for instance, is referred to the S-frame action 
and to the S-frame equations (112)). The reason of such a frame dependence 
is that, in a generic frame, the metric and the dilaton fields are non-trivially 
mixed through the Zr and Z§ coupling functions, so that the associated 
dilaton charge actually controls the matter coupling not to the pure scalar 
part, but to a mixture of scalar and tensor part of the gravi-dilaton field. 

A frame-independent and unambiguous definition of the dilaton coupling 
strengths can be given, however, in the canonically rescalcd Einstein frame 
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(E-frame), where the full kinetic part of the action (114) (including the matter 
and gravi-dilaton sector) is diagonalizcd in terms of the canonically normal- 
ized fields <7 M „, 4> and tpi [13]. Assuming that the dilaton is stabilized by its 
potential, and expanding the Lagrangian term describing the interaction be- 
tween <j> and ipi around the value 4>o which extremizes the potential, we can 
define, in this rescaled frame, the effective masses fhi and charges % for the 
canonical fields ipi. In the weak coupling limit in which Zr ~ Z^ ~ cxp(— <fr) 
one then finds, in particular, that the canonical dilaton charge % deviates 
from the standard "gravitational charge" by the dimensionless factor [13] 



qi ~l 



n i, 



9 lnf4 

z k 



(115) 

4>=<t>o 



For a pure Brans-Dickc model of scalar-tensor gravity one has, for instance, 
q i = 1 (because there is no dilaton coupling to the matter fields in the Jordan 
frame, where dZ % /d<p = 0). For a string model, on the contrary, the coupling 
parameters q { deviate from 1 and are non-universal, in general, since the loop 
form factors Z l tend to be different for different fields ipi. In particular, in the 
conventional scenario which assumes that the loop corrections determining 
the coupling are the same determining also the effective mass of the given 
particle, one obtains large dilaton charges (q i ~ 50) for the confinement- 
generated components of the hadronic masses [48, 49], and smaller charges 
(q { <~ 1) for the leptonic components. In that case, the total dilaton charge of 
a macroscopic body tends to be large (in gravitational units) and composition- 
dependent [55], so that a large dilaton mass (m > 10~ 4 cV) is required to 
avoid conflicting with known gravitational phenomenology. 

This conclusion can be avoided if the loop corrections combine to produce 
a cancellation, in such a way that the value of the coupling parameters q i 
turns out to be highly suppressed with respect to the natural value of order 
one (a scenario of this type has been proposed, for instance, in [56]). In that 
case -C 1, and light dilaton masses (as required, for instance, for a resonant 
interaction with gravitational antennas) may be allowed, without clashing 
with experimental observations. 

In the rest of this section we will focus our attention on this possibility, 
considering the response of the gravitational detectors to a cosmic background 
of massive, non-relativistic dilatons, assuming that the background energy 
density corresponds to large fraction of critical density, and that the dilatons 
are arbitrarily light and very weakly coupled to ordinary matter. 



2.2 Dilaton signals in gravitational antennas 

The operation mechanism of all gravitational antennnas is based on the so- 
called equation of "geodesic deviation" (see e.g. [57]), which governs the re- 
sponse of the detector to the incident radiation. Such an equation is obtained 
by computing the relative acceleration between the world-lines of two nearby 
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test particles, separated by the infinitesimal space- like vector rf , and evolv- 
ing geodesically in the given gravitational background. The interaction with 
a dilaton background can be easily included, in this context, by replacing 
the geodesic paths of the test particles with the world-lines described by Eq. 
(113): one is lead, in this way, to a generalized equation of deviation [55], 

^ + ^Y u V = giJ %v^, (lie) 

which is at the ground of the response of a detector to a background of gravi- 
dilaton radiation (the symbol D denotes covariant differentiation along a curve 
parametrized by the affine time- like variable r). 

This equation implies that a gravitational detector can interact with the 
scalar radiation in two ways: either 

i) directly, through the non-geodesic coupling of its scalar charge to the sec- 
ond derivatives of the scalar background [55, 58]; or 

ii) indirectly, through the geodesic coupling of its gravitational charge to the 
scalar part of the metric fluctuations induced by the dilaton, and contained 
inside the Riemann tensor [59]. 

For a precise discussion of the response of the detector we need to compute 
the "physical strain" h{t) induced by the scalar radiation, which is expressed 
in terms of the so-called "antenna pattern functions" F{6,4>), describing the 
detector sensitivity along the different angular directions. To this purpose, we 
shall rewrite Eq. (116) in the approximation of small displacements £ M around 
the unperturbed path of the text bodies, by setting r\^ = + £ M (r), with 

= const. We then obtain, in the non-relativistic limit, 

C = -L k M k \ (117) 

where 

M k 4 = R k00 1 + qdkd'dp (118) 

is the total (scalar-tensor) stress tensor describing the "tidal" forces due to 
the incident radiation. For the pedagogical purpose of this paper we shall 
assume that the tensor (i.e., gravity- wave) part of the radiation is absent, and 
that the scalar radiation can be simply described as a linear fluctuation of the 
Minkowski metric background and of a constant dilaton background <j>o: 
thus, in the longitudinal gauge, 

ds 2 = (i]^ + 5g^ u ) dx^dx v = (1 + 2^))dt 2 - (1 - 2ip)5 ij dx i dx j , 
^-^o+X, (119) 

so that 

Mij = didjip - Siji> - qdid-jX- (120) 

To discuss the detection of a stochastic background of massive scalar ra- 
diation it is also convenient to expand the fluctuations in Fourier modes of 
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proper momentum p = pn and frequency v = E(p) = (p 2 + to 2 ) 1 / 2 , where 
the unit vector n specifies the propagation direction of the given mode on the 
angular two sphere f2 2 - We obtain 



M y = \ I dp [ d 2 n (2nE) 

^ J -co J Q 2 



TO" 



q^n l n 3 x{p-,n) 



S ij tp(p,n) - n t n 3 ip(p,n) + -j^n^jtpijp,^) 



(note that we are using "unconventional" units in which h = 1, i.e. H = 1/27T, 
for an easier comparison with the experimental variables). We will also assume 
that the dilaton is the only source of scalar metric perturbations, so that 
(p = if) [40]). Introducing the transverse and longitudinal projectors of the 
scalar stresses, defined respectively by 



T- 



5i 



ij TliTlj , 



J ij TliTlj , 



(122) 



defining My = —Fij , and projecting the stress tensor onto the detector tensor 
D y (specifying the geometric configuration and the orientation of the arms 
of the detector), we finally obtain the scalar strain as [58, 60, 61] 



h(t) = D lj F k 



= 5f */ 

z J-oo J n- 



d 2 r 



F sco (n)iP(p,n) 



+ F^(n) X (p,n) 



e 2m(pn-x-Et) _|_ ^ c 



Here 



r- = I T,, • ^/.,, 



P n S — „f— Tji] T 



(123) 

(124) 
(125) 



are the antenna pattern functions corresponding, respectively, to the geodesic 
(or indirect) and non-geodesic (or direct) interaction of the detector with the 
scalar radiation background. 

It should be noted that the scalar radiation, differently from the case of 
the tensor component, contributes to the response of the detector also with 
its longitudinal polarization states. The longitudinal contribution is present 
also in the ultra-relativistic limit to — > 0, p — > E, thanks to the non-geodesic 
coupling (125). In the opposite, non-relativistic limit p — > 0, E — > to, the 
geodesic strain tends to become isotropic, + (m/E) 2 Lij — > Sij, while the 
non-geodesic one becomes sub-leading. 

The results (123) is valid for any type of detector described by the re- 
sponse tensor D*- 7 ', and is formally similar to the expression for the strain 
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obtained in the case of tensor gravitational radiation - modulo the presence 
of different pattern functions, due to the different polarization properties. The 
scalar strain (123) can thus be processed, following the standard procedure, 
to correlate the outputs of two detectors and to extract the so-called signal- 
to-noise ratio (SNR), representing the experimentally relevant variable for the 
detection of a stochastic background of cosmic radiation [62]. 

For our scalar massive background, with spectral energy density 0{p), we 
obtain [58, 60, 61], in particular, 



7 

Jo 



p3 (p2 + m 2 )3 /2 Pi(v / p 2 +to2) p 2 ( % /p2 +m 2 ) 



(126) 



(see also [32] for a detailed computation). Here T is the total (experimental) 
correlation time, N an (irrelevant) normalization factor, Pi and P2 the noise 
power spectra of the two detectors, and j(p) the so-called "overlap reduction 
function", which modulates the correlated signal according to the relative 
orientation and distance of the detectors, located at the positions x 1 and x 2 : 

7(p) = ^ / d2 « Mn) ^(nje 2 "*"-*" 1 -"^. (127) 

The overlap is to be calculated with the geodesic pattern function Ff eo of 
Eq. (124) if we are considering the indirect signal due to a spectrum of scalar 
metric fluctuations, 0^,(p); it is to be calculated with the non-geodesic pattern 
function F™ 9 of Eq. (125) if we are considering, instead, the direct signal due 
to a spectrum of dilaton fluctuations, O x {p). 

We are now in the position of stressing another important difference from 
the case of pure tensor radiation, due to the presence of the mass in the 
noise power spectra Pi. For a typical power spectrum, in fact, the minimum 
level of noise is reached around a rather narrow frequency band outside 
that band the noise rapidly diverges, and the signal (126) tends to zero. As 
v = (p 2 + rn 2 ) 1 / 2 we have, in principle, three possibilities. 

1) If m 3> ^0 then the noise is always outside the sensitivity band Pi(fo), and 
the signal is always negligible. 

2) If m <C v$ then the sensitivity band may only overlap with the relativistic 
sector of the spectrum, for p <~ v <~ 

3) If m ~ i^o, finally, the whole non-relativistic part of the spectrum p <, m 
satisfies the condition Pi{v) ~ Pi(m) ~ Pi{vo). 

It is thus possible to obtain a resonant response to a massive, non-relativistic 
background of scalar particles, provided the mass lies in the band of maximal 
sensitivity of the two detectors [58, 60]. Considering the present, Earth-based 
gravitational antennas, operating between the Hz and the kHz range, it follows 
that the maximal sensitivity is presently in the mass range 



Dilaton cosmology and phenomenology 



39 



1(T 15 eV < m < 1Q- 12 eV. (128) 

Amusingly enough, it turns out that such small values are not so unrealistic 
if the dilaton mass is perturbatively generated by the mechanism of radiative 
corrections. For a scalar particle, gravitationally coupled to fermions of mass 
Mf with dimensionless strength q, there are, in fact, quantum loop corrections 
to the mass of order qMf(A/Mp), where A is the cut-off, which we shall assume 
typically localized at the TeV scale (see for instance [63]). Considering the 
dilaton coupling to ordinary baryonic matter (Mf ~ 1 GeV) the induced 
mass is then: 

™~^V^Wo- 6 eV. (129) 



^lTcv/ V!Gev, 

Thus, a value of q smaller than (but not very far from) the present upper 
limits [53] (imposing q < 10~ 4 in the relevant mass range (128)) is perfectly 
compatible with the possibility of resonant response of the present detectors. 

Quite independently from the possible origin of the dilaton mass, if we 
assume that the mass is in the resonant range (128), and that the bounds on 
q are satisfied, we find that a cosmic background of non-relativistic dilatons 
is possibly detectable by the interferometric antennas of second generation 
- such as Advanced and Enhanced LIGO - provided the background energy 
density is sufficiently close to the saturation of the critical density bound 
[58, 60]. This interesting possibility can be illustrated by considering, for an 
approximate estimate, the simplified situation of two identical detectors with 
Pi = P 2 = P, responding non-geodesically with maximal allowed overlap 
~ g 2 (47r/15) (the numerical factor is referred to the particular case of 
interferometric antennas). Let us suppose, also, that the SNR integral (126) is 
dominated by the peak value Q m of the non-relativistic dilaton spectrum, and 
that such value is reached around p = m (otherwise the response is suppressed 
by the factor (p/m) 4 , [60]). Eq. (126) gives, in this case, 



s^^ q \Tr, (130) 

and the condition of detactable background (SNR > 1) implies 

1 /2 

m 5 / 2p(m) < f h 2 Qm (j^^j x IO- 33 Hz 3 / 2 . (131) 

This condition is compared in Fig. 6 with the expected spectral noise of the 
three LIGO generations (see e.g. [64]), for T = 4x 10 7 s. The region of the plane 
{m, P} corresponding to a detectable background is located above the bold 
noise curves (labelled by LIGO I, LIGO II and LIGO III), and below the dashed 
lines, representing the upper limit (131) for different constant values of the 
parameter q 2 h 2 fi m . This limit may be interpreted cither as a constraint on the 
intensity f2 m , for backgrounds geodesically coupled (q 2 = 1) to the detectors, 
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or as a limit on the non-geodesic coupling strength q 2 , for backgrounds of 
given energy density fl m . As shown in the picture, phenomenologically allowed 
backgrounds are in principle accessible to the sensitivity of next-generation 
interferometers (see also [32] for a more detailed discussion). 
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Fig. 6. The noise power spectra of the three LIGO generations (bold curves), and 
the condition of detectable dilaton background (dashed lines), plotted at different 
values of the parameter q 2 h 2 fl m (ranging from 10" 5 to 10~ n ). 




2.3 Enhanced signals for flat non-relativistic spectra 

The result reported in Eq. (130) is generally valid for a growing spectrum 
with a steep enough slope, as typically obtained in "minimal" models of pre- 
big bang inflation. However, the cross-correlated signal may result strongly 
enhanced with respect to Eq. (130) if the dilaton spectrum is sufficiently flat, 
and if the considered pair of detectors satisfies the condition j(p) — > const 
^ for p 0. 

Let us consider, in fact, the SNR integral (126), which can be written as 

where E = (p 2 + m 2 ) 1 / 2 , and where we can assume that f2(p) is a power-law 
function of p, with an ultraviolet cut-off at p = p±. For a massless spectrum 
(p = E) this integral is always convergent (for any slope), even in the infrared 
limit p = E — > 0: in fact, when p — > 0, the physical strains are produced 
outside the sensitivity band of the detectors, and the noises blow up to infinity, 
Pi{E) — > -Pj(O) — > co. For to 7^ 0, on the contrary, in the infrared limit p — > 
the noises keep frozen at the frequency scale determined by the mass of the 
scalar background, Pi{E) — > Pi(m) = const. In this second case the behavior 
of the integral dependes on j(p) and f2(p). 
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Suppose now that 7(p) — > 70 = const for p — > 0, and that i?(p) ~ p s , for 
p < m. For 5 < 1 we find that the integral is dominated by the infrared limit, 
and gives 



{SNR) 2 ~ 3— / -| ^(p) 

m i Pi{m)P2\m) J p 4 



T ll l, (133) 



m 3 P 1 P 2 



Thus, the integral is infrared divergent [65] for all spectra (even if blue, 6 > 0) 
with S < 1 ! 

This divergence is obviously unphysical, and can be removed by noting 
that the observation time T is finite, and is thus associated to a minimum 
resolvable frequency interval Av = AE = A(p 2 /2m) > T -1 , defining the 
minimum momentum scale 

p min - (2m/T) 1/2 > 0, (134) 

acting as effective infrared cut-off for the integral (133). This implies a modi- 
fied dependence of SNR on the correlation time T in the case of flat enough 
spectra: 

'TV2, 5>1, 



SNR~T^[p*-X_~{ T T l 



5 /\ 6<1. 



(135) 



For S < 1, in particular, there is a faster growth of SNR with T, which may 
produce an important enhancement of the sensitivity to a cosmic background 
of non-relativistic scalar particles, as discussed in [61, 65]. 

It is important to stress that the case j(p) — > 70 = const for p — > has not 
been "invented" ad hoc: it can be implemented, in practice, with detectors 
already existing and operative (or with detectors planned to be working in 
the near future, like resonant spheres). A first simple example, studied in [65], 
refers in fact to spherical, resonant-mass detectors, whose monopole mode 
is characterized by the "trivial" response tensor D % i = 6^. In that case the 
geodesic pattern function (124) is isotropic, 

f ,» = v±^jl, (136) 

pi _|_ m z 

and the geodesic overlap function (127), for two identical spheres with spatial 
separation \xi — x 2 \ = d, is given by 

2 /2p 2 + 3m 2 \ 2 sin(27rpd) 



^ = N\r^ir) - (137) 

This function clearly satisfies the requirement -f(p) — > 70 = const for p — > 0. 

A second example, studied in [61], refers to the so-called "common mode" 
of the interferometric antennas, characterized by the response tensor 
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Df=«V+«V, (138) 

where u l and v % are the unit vectors specifying the spatial orientation of 
the axes of the interferometer. Let us consider, for instance, a geometrical 
configuration where the vectors u and v are coaligned with the X\ and x-i axes 
of a Cartesian frame, respectively, and the direction n ofthc incident radiation 
is specified (with respect to the axes X\, x 2 and x 3 ) by the polar and azimuthal 
angles ip and 9. The computation of the geodesic pattern function (124) gives, 
in that case, 

Ff° = 2- (|) 2 sin 2 0. (139) 

The geodesic overlap function (127), for two coplanar interferometers with 
spatial separation \Ax\ = d, is [61] 




(140) 



where a = 2irpd, and jo, ji, ji are spherical Bessel functions. Thus, also in 
this case, 7 — ► 16n/N = const for p — > 0. 



tT(p) = 



47T 

~N 



3 Late-time cosmology: dilaton dark energy 

In this third lecture we will discuss the possibility that a homogeneous, large- 
scale dilaton field may be the source of the so-called "dark energy" which 
produces the cosmic acceleration first observed at the end of the last century 
[66] , and confirmed by most recent supernovae data [67, 68] . 

Let us recall, to this purpose, that the initial phase of pre-big bang inflation 
is characterized by the monotonic growth of the dilaton and of the string 
coupling g s (see Sect. 1.3): the subsequent epoch of standard evolution thus 
opens up in the strong coupling regime, and should be described by an action 
which includes all relevant loop corrections. Late enough, i.e. at sufficiently 
low curvature scales, the higher-derivative corrections can be neglected, and 
the action can be written in the form of Eq. (114). In that context the loop 
form factors Z((f>), and the dilaton potential V(<j>), may play a crucial role in 
determining the late-time cosmological evolution. 

There are, in principle, two possible alternative scenarios. 

i) The dilaton is stabilized by the potential at a constant value <p = </>o 
which extremizes V{<j>). In this case the loop corrections induce a constant 
renormalization of the effective dilaton couplings (as discussed in Sect. 
2.1), and the Universe may approach a late-time configuration dominated 
by the dilaton potential, with H 2 ~ V((f> a ). 
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ii) The dilaton fails to be trapped in a minimum of the potential, and keeps 
running even during the post-big bang evolution. In this case the late-time 
cosmological evolution is crucially dependent on the asymptotic behavior 
of the factors Z{<f>). 

These two different possibilities have different impact on the so-called "coin- 
cidence problem" (i.e. on the problem of explaining why the dark-matter and 
dark-energy densities are of the same order just at the present epoch), as we 
shall discuss in the following subsections. 

3.1 Frozen dilaton in the moderate coupling regime 

The first type of scenario can be easily implemented [69] using a generic non- 
perturbative potential which is instantonically suppressed (V ~ exp(— 
in the weak coupling limit g 2 — > 0, and which develops a non-trivial structure 
with a (semi-perturbative) minimum g 2 ~ olgut ~ (M s /Mp) ~ 0.1 — 0.01 in 
the regime of moderate string coupling. A typical example is the "minimal" 
potential given, in the E-frame, by [70] 



V{4>) = m 2 v e* 1 ^-* 1 ) + /fcr* 2 ^ - * 1 ) 



= exp[-7(0-0i)] 



(141) 



where k\, k 2 , /3, e, 7 are dimensionless parameters of order one (see Fig. 7). 
The presence of a local minimum at cf> ~ 4>i allows solutions with <f> = const 
during the radiation-dominated phase, and (for appropriate values of my) 
may also lead to a late phase of accelerated expansion driven by the potential 
energy V((j>o), provided the dilaton is not permanently shifted away from the 
minimum 4>q by the transition to the matter-dominated epoch [69] . 




Fig. 7. Plot of the potential (141) for ki = k 2 = /3 = 7 = 1, e = 0.1, cj>i = -3, 
my ~ 0.1, and a local minimum (independent of my) at <j>o = —3.112, corresponding 
to g 2 s = exp(0 o ) ^ 0.045. 
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Let us consider, in fact, the equation of motion of a homogeneous dila- 
ton field <p(t) in the conformally rescaled E-frame (with metric g) where the 
graviton kinetic energy is canonically normalized, and let us assume that the 
rescaled matter sources can be described as a perfect fluid of energy density 
p, pressure p, and dilaton charge a. Starting from an action of the type (114) 
we find that the generalized dilaton equation, for a cosmological background, 
takes the form 

A(<f>) + 3H<p) + B{4>)tf + ^ + A| [C(<f>) (p-3p)+a} = 0, (142) 

where A, B and C are functions describing the rescaled (E-frame) loop correc- 
tions. For a minimally coupled field, for instance, A = 1, B = C = a = 0; for 
the dilaton, at tree- level in the string coupling, A = C = 1, B = 0. In the most 
general case we find that a stable dilaton configuration with <fi = = <j> is pos- 
sible, in the radiation era (p = 3p), if the scalar charge of the fluid is negligible, 
a = 0, and the dilaton extremizes the E-frame potential, dV/d(j) = 0. 

When the Universe becomes matter-dominated (p = 0), however, a new 
acceleration <fi = —A~ 1 \pCp is suddenly generated, which tends to remove 
the dilaton away from its equilibrium position. Such an acceleration is in 
competition with the restoring force cf) = — A~ l {dV /d<j>) (see Eq. (142)). The 
possibility that the dilaton may bounce back to the stable minimum cf> — </> , 
driving the Universe towards a final phase of accelerated, potential-dominated 
expansion, thus crucially depends on the values of two parameters: the (loop- 
corrected) strength \pC((f>o) of the dilaton coupling to dark matter, and the 
slope of the dilaton potential (141), determined by the mass scale ray which 
also controls the amplitude of the minimum, V(<fto) ~ m v - Such an amplitude, 
on the other hand, should correspond to the present Hubble scale (U(</>o) ~ 
Hq), in a realistic model able to describe the present phase of accelerated 
expansion. 

It can be shown, with a simple numerical analysis, that the values of 
the coupling strength allowed by present gravitational phenomenology are 
compatible with a late-time phase dominated by the potential only for a finite 
range of values of V((po), depending on the value of the dilaton coupling at 
the equality epoch [69]. Using the phenomenological upper limit |C cq | — 0.1 
one finds that the dilaton, after a smal shift at t = t cq , bounces back to the 
minimum provided 10~ 7 _ff O q ^5 itiv ^ -^cq (which includes the realistic case 
my ~ H ~ 10~ 6 H eq ) [69]. Smaller values of |C oq | correspond to a larger 
mass interval. We can say, therefore, that the coincidence problem (i.e., why 
V(ipo) ~ Hq), in this context remains, but is somewhat alleviated because 
- thanks to the dynamical correlation between the amplitude V{<j>o) and the 
matter-dilaton coupling - only a restricted range of values is allowed for U(</>o). 
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3.2 Running dilaton: saturation of the loop corrections and 
asymptotic "freezing" 

The second possibility, which will be discussed here in more detail, in the case 
in which the dilaton is not stopped by the structures formed by the potential 
around g 2 = 1, and keeps rolling towards +oo along a smoothly decreasing 
potential. A possible example of non-perturbative potential of this type is 
given, in the E-frame, by [71] 

V = c{m 2 v ( ^ \ \ [e- ft oxp (-« - e- ft cx p<-« 1 , (143) 

where b\, c±, 0i, fti are dimensionless parameters, with < f3\ < 02. This 
potential is instantonically suppressed in the weak coupling limit <j) ~ > ~ oo, 
and is exponentially decaying as 

V = m 2 v {p 2 - )9i)e-* + G(e- 2 ^) (144) 

in the limit <f> — * +oo (see Fig. 8). In this case, as we shall see, we can obtain a 
scenario of "coupled quintessence" [72] in which the late Universe approaches 
a (possibly accelerated) state dominated by a mixture of kinetic and potential 
energy density, and the coincidence problem may find a satisfactory solution 
thanks to the dilaton-dark matter interactions. 
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Fig. 8. Plot of the potential (143) for bi = 1, ci = 10, f3i = 0.1, jh = 0.2 and 
my = 1. The dilaton is monotonically growing from the string perturbative vacuum 
along a "bell-like" non-perturbative potential. 



In this case, however, a realistic scenario requires some mechanism of sat- 
uration of the loop corrections, so as to keep the present effective values of 
gravitational and gauge couplings approximately constant and sufficiently 
"small", even in the large "bare coupling" limit 4> — * +°°- As discussed in 
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[73], such a saturation can be obtained thanks to the large number of fields 
(e.g. gauge bosons) entering the loop corrections, assuming (as in models of 
"induced gravity") that the loop form factors of Eq. (114) have a finite limit 
for (j) — > +oo, and that can be approximated by a Taylor expansion in powers 
of the inverse bare coupling g 2 = exp <f>. Applying these assumptions to the 
gravi-dilaton form factors, to the potential, and to the dimensionless param- 
eters qi{4>) controlling the dilaton charge-density of the various matter fields, 
we can set, for cf) — ► +oo, 

ZrW = c? + &ie"* + 0(e- 2 ^), 

q i {4>)=q Qi +0{e- 2 +). (145) 

The dimensionless coefficients c 2 and c 2 , of this expansion are typically 
of order N <~ 10 -2 , because of their quantum-loop origin and of the large 
number N of gauge bosons in GUT groups like E$. This is in agreement with 
the fact that c\ controls (according to the action (114)) the asymptotic value 
of the ratio between the string and the Planck length scale, c\ = (A s /Ap) 2 , 
which is indeed expected to be a number of the above order. The coefficients 
&i, &2 ■ • - j on the contrary, are numbers of order one. Note that the expansion 
of V{4>) agrees with the asymptotic form of the potential (144). 

We should note, finally, that the asymptotic values of the dilaton charges, 
qoi, have to be strongly suppressed for the ordinary components of matter 
(such as baryons) and for electromagnetic radiation: if we want a dilaton field 
active on a cosmological scale of distances, in fact, we need long-range inter- 
actions, and we must avoid unacceptable deviations from the standard gravi- 
tational phenomenology by suppressing the dilaton couplings, as discussed in 
Sect. 2.1. For the (possibly exotic) components of dark matter, however, there 
is no strict phenomenological bound imposing such suppression: in that case, 
the asymptotic charge qo could be nonvanishing, and of order one, leading to 
interesting late-time deviations from the standard cosmological scenario. 

For a simpler illustration of this possibility it is convenient to work in the 
diagonalized E-frame, obtained from the metric g of Eq. (114) through the 
rescaling 

g l w = c l Z R 1 9»v ( 146 ) 
The action (114) becomes, in this new frame 



-n+h 2 ^) (W) 2 -V{<t>) 



S m (g, (j), matter), 

(147) 



where 

2 



d\nZ R \ ~ 4 2 . 



* W= 3 (-«T^) - 2 ^. V{<t>)=c\Z^V. (148) 
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Assuming that the matter action S m describes a perfect fluid with a dark- 
matter component p m , a baryon component pi,, and a radiation component 
p r = 3p r , the cosmological Einstein equations for the action (147) can then 
be written (omitting the tilde, and in units 2A| = 1) as: 

6H 2 = p r + Pb+ Prn + P<j>, 

4H + 6H 2 = -^-- P4> , (149) 

where 

k 2 U) ., k 2 (cb) ., , 

P* = -^V + V, p* = -^V - V. (150) 

The associated dilaton equation, assuming a negligible density of dilaton 
charge for baryons and radiation (ay = = cr&), can be written as [71] 

k 2 + 3Hj>) + kk'tf + V + X - W ( Pb + Pm ) + o m \ = 0, (151) 

where we have defined ip = —ItlZr, and the prime denotes differentiation 
with respect to <p. The combination of Eqs. (149)— (151) leads, finally, to the 
equations of energy- momentum conservation for the various fluid components: 

p r + AHp r = 0, 

ib' ■ 

p b + 'iHp b - Y^P b = ' 

p m + 3Hp m ~Y^Pm- -y-0 = 0, 

P4, + 3H( P<t) + P<p ) + ^<f> [th' { Pb + Pm ) + a m ] = (152) 

(the last equation is simply the dilaton equation (151), rewritten in fluido- 
dynamical form). 

Let us now concentrate on the coupled dark-matter/dilaton system, and 
note that there are two types of interactions between these two cosmic sources: 
a first one, specific to the particular type of dark matter field, generated 
by the "intrinsic" dilaton charge a m ] and a second one, more "universal", 
generated by the standard dilaton coupling to the trace of the stress tensor, 
and associated to the tp' terms of the above equations. Both types of coupling 
are renormalized by the loop corrections, but with opposite effect according 
to the asymptotic limits of Eq. (145). In fact, the dilaton charge tends to 
grow, and to reach a constant asymptotic value as </> — > +oo. The coupling 
parameter ip' , on the contrary, tends to be exponentially suppressed as 

^' = -(lnZ fl )'->^A 0^+oo. (153) 



As a consequence, after the transition to the matter-dominated phase, the 
Universe may enter two different types of dynamical regimes [71]. 
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1) If the dark-matter charge a, m is still negligible at the beginning of the 
matter-dominated phase (as well as the dilaton potential, expected to become 
important only near the present epoch), then the Universe enters the so-called 
"dragging regime", in which p m is coupled to <fi through the ip' terms of Eqs. 
(153), and the evolution of the (still subdominant) dilaton kinetic energy p$ 
is "dragged" by p m . 

The cosmic evolution, during this regime, can be analytically described (in 
an approximate way) by noting that the loop factor k(<p) goes to a constant 
at late enough time scales, 

k((f>) -» k = V2 — , -» +oo, (154) 

Cl 

according to Eqs. (148) and (145). Introducing the canonical variable (f> = k (f) 
(see the action (147)), and neglecting the subdominant contributions of p r and 
Phi we can then rewrite the coupled equations (151), (152), for the dragging 
regime, as follows: 

+ 3^0+1^ = 0, (155) 
p m + 3Hp m - ^p m (j> = 0, (156) 

where e = ip'/ko ~ e _ ^/(\/2ciC2) < 1 is the effective coupling parameter. 

Neglecting the time dependence of e with respect to that of H and 4> (for 
small enough time intervals), we find that the system of equations (149), 
(155), is satisfied by 

$~-2eH. (157) 

Thus, from Eq. (156), 

• 2 

Pm ~ a~ (3+e2) ~H 2 ~j ~p^, 

a~t 2 /(3+« a ). (158) 

During this phase the dark-matter and the (kinetic) dilaton dark-energy densi- 
ties are characterized by the same time-evolution, which slightly deviates from 
the standard behavior of a dust-dominated Universe (p ~ a~ 3 , a <~ i 2 / 3 ). The 
kinematics, however, remains decelerated (as e«l). 

2) A second, possibly accelerated, "freezing regime" is eventually reached 
in the limit in which the dilaton potential comes into play, and the coupling 
induced by the intrinsic charge density er m becomes dominant with respect to 
the exponentially suppressed coupling due to ip 1 . 

Using again the canonical variable 0, assuming that a m = q{4>)p m (for 
a homogeneous fluid), and considering the asymptotic limits q((f>) — ► q , 
V = Vq cxp(— (j>) of Eq. (145), we can rewrite the coupled dilaton-dark matter 
equations (152), for the freezing regime, as follows: 
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p m + 3Hp m - -Y-pm4> = 0, 

p^ + 6H Pk + ^- Pm = O. (159) 

We have defined the kinetic and potential energy densities, p k and pv, respec- 
tively as 

Pfe = y, p v = V(j) = V e-^ ko , p* = p k +Pv. (160) 

The system of equations (159), (149) (with p r = pb = 0) can be solved by a 
late-time configuration in which p m , p<t„ V and H 2 scale in time in the same 
way, so that the critical fractions of dark-matter density, Q m = p m /6H 2 , 
dilaton kinetic energy, Q k = p k /6H 2 , and potential energy, ft v = V/6H 2 , are 
separately frozen at constant values determined by k and q only (i.e. by the 
parameters ci, c 2 , qo of the asymptotic expansion (145)). A simple analysis 
gives [71] 

3k 2 _ 3k 2 + g (g + 2) 

4 " (9o + 2)2' "V- {qo + 2) 2 > 

ftcf, = ^k + &m = 1 — &4>, (161) 

where feo is given by Eq. (154) (see also [32] for a detailed computation). 

In this asymptotic state the Universe is thus dominated by a fixed mixture 
of dark-matter and dilaton (kinetic plus potential) energy density. The dilaton 
fluid has equation of state 

= P<p = - Oy = qojqo + 2) 
W P<p n k + Qv 6kl + qo{qo + iy 

and can play the role of the dark energy fluid responsible for the observed 
cosmic acceleration, provided q > 1. 

In fact, by rewriting the Einstein equations (149) for H in the form 

2H 

l + T^ = f2v-O k , (163) 

we obtain 

d H 3 , „ „ ^ 1 <7n — 1 , 

^ = 1 + ^=2^-^-2 = ^- < 164 > 

The expansion is accelerated (a > 0) for q > 1 or q a < —2. The second case 
(corresponding to an acceleration of superinflationary type, with H > 0) is to 
be excluded, however, in our context, as it would imply fi m < according to 
Eqs. (161). Thus, acceleration is only possible for q > 1. The explicit form 
of this asymptotic solution can be finally obtained through the integration of 
Eq. (164), which gives 
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a ~ i(*+2)/3 



H ~ a 



-3/(90+2) 



(165) 



from which 



■2 



Pm 




Vbe 



-<t>/ko 



-6/( g o+2) 



(166) 



To illustrate the smooth background evolution from the initial radiation 
phase to the intermediate dragging phase, and to the final freezing regime, we 
shall conclude this subsection by presenting the results of an exact numerical 
integration of the string cosmology equations (149)-(152). For our illustrative 
purpose we will assume that Zr and Z^ are given by the expansion (145) 
truncated to first order in exp(— <j>), with b\ = b<i = 1, cf = 100 and c 2 . = 30. 
We will adopt the model of dilaton coupling already used in [71] , parametrized 
by the time-dependent charge 



with c 2 = 150 and q n = 2.5. We will also use the E-frame potential (144), 
with /?i = 0.1, 02 = 0.2, and c\my = 10~ 3 iJ cq . The last choice, which implies 
my ~ Ho, is crucial to obtain a realistic scenario in which the asymptotic 
accelerated regime starts at a phenomenologically acceptable epoch (see [71, 
32] for a discussion of the mass scale of the non-perturbative dilaton potential, 
and of the degree of fine-tuning possibly required for realistic cosmological 
applications). Finally, we will integrate our equations imposing the initial 
conditions p^(ti) = p r (ti), p m {ti) = 10~ 2< Vr(*i) 5 p b (ti) = 7 x 10~ 21 p r (ti), 
tpiti) = -2, at the initial scale H{U) = 10 40 iJ cq . 

The obtained scaling evolution p = p(a) is illustrated in Fig. 9 for the 
various cosmic components. We can note that, at large enough times, baryons 
(full line) and radiation (dotted line) arc fully decoupled from the dilaton, and 
obey the standard scaling behavior (p r ~ a , p b ~ a~ 3 ). The late-time dark- 
matter evolution, on the contrary, is closely tied to the dilaton evolution, 
and the ratio of their energy densities becomes asymptotically frozen at a 
constant. With the particular numerical values used in this example we obtain 
an asymptotic configuration characterized by ~ 0.73 and Q m ~ 0.27, with 
a dark-energy equation of state w ~ —0.76. 

3.3 Non-local coupling and pressure backreaction 

Another interesting asymptotic configuration can be obtained in the case in 
which the dilaton is non-locally coupled to the dark matter components, as 
discussed in subsection 1.2. In that case, the fractions fl m and fi^ may also 
become frozen at constant asymptotic values, but the background evolution 
turns out to be decelerated for any values of q , as the dark matter develops 
an effective pressure which tends to compensate the accelerating action of the 



q{4>) = 5o 



c' 



(167) 



c 2 _|_ e 9o0' 
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Fig. 9. Late-time evolution (on a logarithmic scale) of the various components of 
the cosmic energy density. The plots are the result of a numerical integration of Eqs. 
(149)-(152). 

dilaton potential. This effect is new, and will be illustrated in some details in 
this subsection. 

We start assuming that, in the matter part of the action (114), the dilaton 
is non-locally coupled to the sources through the variable £(</>), as in the action 
(28). However, differently from the action (28), we will assume (for simplicity) 
that the dilaton potential is local, V — V(<fi): the gravi-dilaton part of the 
action is thus identical to that of Eq. (114), and our model is described by 



Let us vary this action with respect to g and 0, and evaluate the resulting 
(general covariant) field equations in the limit of a homogeneous, isotropic, 
spatially fiat background, using the results of Subsection 1.2. We obtain a set 
of equations similar to Eqs. (45)-(47) for what concerns the dilaton charge 
density a(4>), but different for the potential (which now is local), and for the 
presence of the loop corrections Zr, Zq. 

Let us finally transform the equations in the E-frame (using the rescaling 
(146)), and consider the asymptotic limit in which p r , pb are negligible, and 
the dark matter is coupled to the dilaton only through its intrinsic dilaton 
charge (namely, the limit in which ip' ~ 0) . The resulting equations (omitting 
the tilde, in units 2Ap = 1) can be written as 




(168) 



6H 2 = 



Prn + P4>, 



(169) 
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4H + 6H 2 = -p -^p, (170) 

p m + 3if ( Pm + ^) - ^0 = 0, (171) 

p + 3H( /90 +p ) + ^ = O, (172) 

with andp«/> defined by Eq. (150), as before. A comparison with the asymp- 
totic limit of Eqs. (149), (152), shows that the genuinely new effect of the 
non-local interactions is the appearance of an effective pressure term a m /2 
for the dark matter component. Indeed, the new terms present in Eqs. (170), 
(171), can also be obtained from the standard Einstein equations through the 
shift p m — > p m +a m /2. 

We are now in the position of asking whether or not this modification 
(of non-local origin) may change the results of the previous subsection, in 
particular those concerning the asymptotic freezing configuration. We shall 
consider, to this purpose, the limit in which a m — > qop m and V = Vo exp(— (p), 
using the canonical variable <f> as in the previous computations. 

Let us look for solutions of Eqs. (169)-(172) by requiring for p m , pk and 
pv the same scaling behavior, and thus imposing, as a first condition, that 

EHL = P±. (173) 

Pm P4> 

Using Eqs. (171), (172) for p m and p^, and the Einstein equation (169), we 
obtain 

$_ _ 6fc 

We are denoting with a bar the fractions of critical energies for the new 
configuration associated to the non-local equations, to distinguish it from the 
"local" freezing solution of Eq. (161). We also impose, as a second condition, 
that 

ElIL = Py. (175) 

Pm PV 

The definition (160) of pv 1 together with Eq. (171), gives then 



jj = 3fc , (176) 
which, combined with eq. (174), leads to 

2 + qo 9o + 2 

From the definition of and Eq. (176), on the other hand, we have 
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The insertion of this result into Eq. (177) finally gives 

= 3*3(2 - «,) +4«, 

4(2 + g ) V ' 

The combination of Qk and fly provides now the values of fl^, Q m , and 
the equation of state w, according to the definitions (161), (162). As we are 
interested in the kinematical properties of the solution we shall compute, in 
particular, the acceleration parameter aj (aH 2 ): dividing by 6H 2 the modified 
equation (170) we obtain 

jL = l(j2 v -72 k )-p T2 m -^ (180) 



from which 



! + 4 = 4. ( 181 ) 



aH 2 ~ H 2 2' 



quite independently of the values of k and q l The integration of H finally 
provides a <~ t 2 / 3, and H 2 ~ p ~ a~ 3 , as in the standard phase of dark-matter 
dominated evolution. 

The considered model of non-local coupling is thus associated to an asymp- 
totic freezing phase which is decelerated, and in which the dilaton energy den- 
sity has the same dynamical behavior of a dust fluid, p$ ~ p m <~ a~ 3 , in spite 
of a pressure which is nonvanishing, in general: 

3/bn 



w = % o 7° - . (182) 
2 3k 2 + q 



This result can be understood by noting Eqs. (180) and (181) together imply 

77 fc -77 v + ^77 m = _L(p, + ^p) =o, (183) 

namely a zero total pressure for the coupled dilaton-dark matter fluid (see Eq. 
(170). The dark matter pressure associated to the non-local effects thus gen- 
erates a backreaction which exactly compensates - at least in this model - the 
dilaton pressure, leading the system to restore, asymptotically, the standard 
dust matter configuration. 



3.4 Main differences from uncoupled models 

Let us come back to the class of models in which the dilaton is locally coupled 
to the dark matter components, as discussed in subsection 3.2. If we identify 
the accelerated freezing phase with our present cosmological phase, and thus 
the energy density of the dilaton field with the "dark-energy" density respon- 
sible for the present cosmic acceleration, we are lead to a dilaton model of 
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dark energy which is substantially different from the conventional models of 
quintessence [74] based on a rolling scalar field, uncoupled to dark matter. 

A first, important (conceptual) difference concerns the mentioned prob- 
lem of the cosmic coincidence. In the considered class of dilaton models this 
problem, if not solved, is at least relaxed: in fact, the dark-energy and dark- 
matter densities are of the same order not only today but also in the future 
(forever), and also in the past for a significantly amount of time, depending 
on the beginning of the freezing epoch (see below) . 

A second, more phenomenological difference concerns the scaling behavior 
of the baryonic and dark-matter components of the dust fluid during the 
freezing epoch. Because of the coupling to the dilaton, the dilution in time 
of the dark-matter density p m is slower than the standard baryon dilution, 
pb ~ a~ 3 : in particular, the ratio pb/p m decreases in time as 

£L „ a -3W(2+9o) Q84) 
Pm 

(see Eq. (166)). This could explain why the present fraction of baryons is 
small (~ 10~ 2 ) in critical units - provided the accelerated epoch has an early 
enough beginning. Direct /indirect measurements of the past value of the ratio 
Pb/Pm, compared with its present value, could provide unambiguous tests of 
this class of models. 

Finally, concerning the beginning of the accelerated epoch, it is important 
to stress that in dilaton models the acceleration can start much earlier than 
in models of uncoupled quintessence [75, 76]. 

For a simple illustration of this point we may consider a model in which, 
during the accelerated regime, there are two types of sources with different 
dynamical behavior: i) an uncoupled component p u , with pressure p u = and 
scaling behavior p u <~ a -3 (represented by baryons and, possibly, by a fraction 
of non-baryonic dark matter uncoupled to the dilaton); ii) a coupled compo- 
nent p c , with pressure p c — wp c , and a slower scaling behavior p c ~ a~ 3 ( 1+tu ) 
(represented by the dilaton and by the fraction of dark matter coupled to the 
dilaton). Thus, even if today p c dominates, and drives an accelerated evolu- 
tion, at ealy enough times the Universe was dominated by p u , and decelerated. 
From the Einstein equations 

6H 2 = p u + p c , 

AH + 6H 2 = -p c = -wp c , (185) 
we obtain that the acceleration switches off at the scale a acc such that 

= 4[^-(i + 3»)(fl B -i)] acc = o, 



acc 




aH 2 V \ H 2 2 



(186) 

where fi u = p u /QH 2 , fi c = 1 — Q u . In terms of the present values fi® of 
these fractions the above condition becomes 
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/ \ — 3 / \ — 3(l+ui) 

— ) = (l + 3«,)(/2°-l)(^) , (187) 

\ a o / \ a o / 

and fixes the beginning of the acceleration at the redshift scale z acc such that 

-l/Sw 

-1. (188) 

Consider now a model of uncoupled quintessence, in which the uncoupled 
component corresponds to the totality of the dark-matter fluid (plus subdom- 
inant contributions), i.e. L?2 = L?J^. Using the recent analysis of the SNLS 
Collaboration [68] , based on present observations of supernovae la and large- 
scale structure, one finds 0.2 < < 0.4, and —1.2 < w < 0.8. One then 
obtains, from Eq. (188), 0.4 < z acc < 1 (see Fig. 10, left panel). 

If we consider instead a model of dilaton dark energy, then the uncoupled 
component may range from the baryon component fib to some fraction of 
the dark matter component fl m . In the "maximally coupled" version of the 
model, in which f2^ = {]®, we can re-apply the supernovae results of SNLS 
with f?° ~ 0.04 - 0.05, to obtain w ~ -0.65 ± 0.15. Eq. (188) then implies 
[76] z acc ~ 3 — 4 (see Fig. 10, right panel). 



no 

&acc 



(1 + 3w) 




-1.4 -1.2 -1 -0.8 -0.6 -0.4 -1 -0.9-0.8-0.7-0.6-0.5-0.4 



w w 

Fig. 10. Beginning of the accelerated epoch for dark-nergy models with uncoupled 
(left) and fully coupled (right) dark matter, according to the observations of type la 
supernovae in the SNLS dataset. The plotted curves are obtained from Eq. (188), 
for constant values of the present fraction of the uncoupled dust matter fi°. 



Thus, dilaton models of dark energy are compatible with a beginning of 
the cosmic accelerations at epochs much earlier than those suggested by other 
models, according to the most recent supernovae data. The extension back in 
time of the accelerated regime might have a significant impact on the dilution 
of baryons, according to Eq. (184). Finally, strongly coupled models tend to be 
compatible with a "less negative" parameter w (see Fig. 10), thus alleviating 
the need for "phantom" dark energy [77] with "supernegative" (w < — 1) 
equation of state. 
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Appendix A 

In this Appendix we present a detailed derivation of the equations of mo- 
tion (29), (34), starting from the action (28) which includes non-local dilaton 
interactions. 

The functional derivation of the action with respect to the metric g' ll '{x) 
contains, besides the standard contributions leading to Eq. (3), the new non- 
local contributions V^ v {x) and L flv (x), and can be written as follows, 



(V and C' m are defined by Eq. (33)). We need now to compute the functional 
derivation of exp(— £). Using the definition (25) we obtain 




(A.l) 



J X 



where 




(A.2) 



(A.3) 




(A.4) 



where 7 M „ is defined in Eq. (30). Thus: 
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V IMf = 



.1 



2X a s 



(A.5) 
(A.6) 



where Iy and I m are defined in Eqs. (31), (32). Inserting these results in Eq. 
(A.l), multiplying by (— 2Af _1 ) exp(—^))/^—g, and imposing the condition of 
zero functional derivative, one is finally lead to Eq. (29). 

Let us now consider the functional derivative with respect to <fi(x). Sepa- 
rating the local and non-local terms, as before, we obtain 



5<P(x) 2\t x 

+ A(x) + B(x), 



[R + 2V 2 q)-(Vq)) 2 + V] i 



where 



B(l ) = /d^V(^C)„^. 



-«(*') 



(A.7) 

(A.8) 
(A.9) 



The functional derivative of the variable (25) leads to 



-«*') - 



8(j){x) 



+ (V^g^V^W) # (fa - <j) y ) [8 d +\x - x 1 ) - 8 d+1 (x - y) 



-d u 



8{4>x> ~ 4>y) 



S d+1 (x-y) 



(A.10) 



where = d/dy^ 7 and 8' denotes the derivative of the delta function with 
respect to its argument. 

The first term of this integral exactly cancels the term containing d^e~^ 
in the last part of the integral; also, the third term exactly cancels the term 
containing d fl [8((j) x > — <p y )\ in the last part of the integral. Thus, we are left 
with: 

<(x>) = 8 d +\x-x>) J dd+ly ( V z^-* v £(v^^ - <p y ) 



8<p(x) 



-gd^ 



(A.ll) 
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The second term on the right-hand side of the above equation can be conve- 
niently rewritten as 



,-<t> 



-(J 



(A.12) 



For the first term containing 5' we can exploit the properties of the delta 
function, and the identities 



dy = 7-d^, 

<Py 



d 1 d 

d<j) y $ y dy ' 



(A.13) 



to obtain 



K d ( d d y (V~9^V^W) 6' fa - 4>y) 



= a: 



= -e~« (a;) +e^ (a;) J(x), (A.14) 

where J is the integral defined in Eq. (35). Inserting the results (A.12), (A.14) 
in Eq. (A. 11), using the definitions of A and B, and integrating over x', we 
finally obtain 



A{x) + B(x) = 
+ e 



-5e 



2A s d " 1 



gC' m ) ( e -«- e -*J) a 



2A S ° 



T[ Iv — In 



(A.15) 



Summing this contribution to the other terms of Eq. (7), multiplying by 
(2Af _1 exp <p/\/—g), and imposing the vanishing of the functional derivative, 
we are lead to the equation of motion (34) for the dilaton. 
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